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Preface 


I  selected  my  thesis  topic  for  three  reasons.  The  primary  reason  was  to  apply 
mv  newly  acquired  classroom  knowledge  to  a  practical  fluids  problem.  Secondly, 
I  wanted  to  do  research  in  a  field  that  could  benefit  current  Air  Force  interests, 
and  finally,  to  do  an  experimental  thesis  to  integrate  thirteen  years  of  electronics 
experience  with  my  newly  acquired  aeronautical  fluid  dynamics  knowledge. 

The  successful  completion  of  any  project  requires  the  hard  work  of  many  peo¬ 
ple.  Of  all  the  hard  working  people,  Dr.  W.  C.  Elrod,  my  thesis  advisor,  played 
the  most  important  part.  I  am  deeply  grateful  for  his  guidance  and  support.  He 
was  always  available  and  many  hours  were  spent  in  guiding  me  to  a  successful  con¬ 
clusion  to  my  research.  I’d  also  like  to  thank  my  committee  members,  Dr.  M.  E. 
Franke  and  Lt  Col  P.  I.  King,  for  the  benefit  of  their  extensive  knowledge  and  helpful 
suggestions.  Their  support  was  greatly  appreciated. 

Someone  once  told  me  a  good  sponsor  is  one  that  identifies  the  general  scope 
of  the  project,  provides  the  necessary  funds,  and  then  leaves  the  researcher  to  do  his 
job.  If  that’s  so,  then  I  was  truly  blessed!  My  thanks  go  out  to  the  AFWL  team, 
especially  Capt  Gregory  Gates,  for  suggesting  the  project,  providing  the  necessary 
information,  and  last  but  not  least,  providing  the  funding  to  put  it  all  together.  This 
report  is  just  a  small  piece  of  the  big  picture,  but  with  continued  support,  it  could 
also  be  the  first  piece  in  the  solution  to  a  bigger  puzzle. 

Laboratory  work  requires  the  sucessful  acquisition  and  operation  of  many 
pieces  of  hardware.  Mr.  Nicholas  Yardich  and  the  AFIT/ENY  technicians  were 
always  readily  available  with  their  knowledge  and  guidance.  I  especially  thank  Mr. 
Yardich  for  all  the  help  in  ordering  equipment,  Mark  Derriso  for  always  finding  the 
part  I  needed,  and  Jay  Anderson  for  knowing  how  to  put  it  all  together  and  make  it 
run. 
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Experimental  work  requires  the  design  and  fabrication  of  many  pieces  of  hard¬ 
ware,  and  the  AFIT  Model  Shop  provides  the  expertise  necessary  to  accomplish  the 
job.  At  one  time  or  another,  every  one  of  the  craftsmen  at  the  model  shop  played  a 
role  in  the  fabrication  of  hardware  for  my  research  and  the  high  quality  of  their  work 
is  greatly  appreciated.  I  want  to  especially  thank  Jack  Tiffany  for  knowing  where 
to  find  the  difficult-to-obtain  parts  and  John  Brohas  for  doing  such  an  excellent  job 
on  the  test  section  and  diffuser.  John  Brohas’s  technical  experience  combined  with 
his  craftsman  skill  were  significant  factors  in  establishing  and  fabricating  the  final 
design  in  a  timely  manner. 

Experimental  work  is  taking  theory  and  putting  it  to  work.  During  the  appli¬ 
cation  of  the  theory,  the  experimentalist  identifies  problem  areas  and  attempts  to 
determine  solutions  to  the  problems.  Many  of  the  problems  have  been  seen  before 
but,  unless  an  experimentalist  has  had  past  experience,  he  must  sometime  reinvent 
the  wheel.  Capt  F.  Tanis’s  suggestions  and  gentle  guidance  based  on  his  experimen¬ 
tal  experience  and  computer  knowledge  saved  me  a  tremendous  amount  of  time  and 
was  greatly  appreciated. 

I  would  also  like  to  thank  my  wife  and  three  children  for  their  tremendous 
patience  and  support  throughout  this  entire  program.  I  would  not  be  where  I  am 
today  if  it  wasn’t  for  their  encouragement. 

Finally,  I  would  like  to  dedicate  this  paper  to  my  father  and  mother,  Stephen 
W.  Stiglich,  Sr.  and  Norma  R.  Stiglich.  Their  continued  belief  in  me  has  never 
faltered.  Thanks  Mom  and  Dad! 


Stephen  Walter  Stiglich,  Jr. 
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Abstract 

A 

Chemical  lasers  require  a  cavity  that  establishes  and  maintains  the  proper  gas 
dynamic  properties  during  lasing.  The  design  and  performance  of  a  flow  system 
capable  of  supporting  the  hypersonic  flow  conditions  in  a  lasing  cavity  are  described. 
Using  cold  air  as  the  working  medium,  the  flow  control  system  configuration  and 
nozzle-cavity-supersonic  diffuser  assembly  configuration  were  developed  to  establish 
acceptable  flow  conditions  in  the  test  section.  Performance  evaluation  was  based 
on  pressure  measurements  in  the  nozzle-cavity-diffuser  assembly  and  schlieren  pho¬ 
tographs  of  the  flowfield  in  the  cavity. 

Flow  conditions  in  the  test  section  were  broken  up  into  three  different  regions: 
flow  in  the  hypersonic  nozzles,  flow  in  the  base  region  and  flow  in  the  cavity  region. 

Flow  in  the  nozzles  was  analyized  using  one-dimensional,  steady,  isentropic 
flow  theory.  Test  results  indicated  that  the  hypersonic  nozzles  performed  to  design 
specifications. 

The  Korst  two-dimensional  base-pressure  flow  model  was  used  to  describe 
the  flow  in  the  nozzle  exit  plane  and  base  region.  Experimentally  calculated  Mach 
numbers  and  static  pressures  corresponded  very  closely  to  theoretical  values. 

Static  pressure  ports  and  schlieren  photographs  were  used  to  describe  the  flow- 
field  conditions  in  the  cavity  region.  Pressure  measurements  indicated  that  super¬ 
sonic  conditions  were  reached  in  the  cavity  for  specific  supersonic  diffuser  throat 
areas  settings,  but  conditions  were  short  lived.  Boundary  layer,  frictional,  and 
three-dimensional  effects  were  suspected  as  the  main  contributors  to  the  flowfield 
degradation.  -  f 
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EXPERIMENTAL  INVESTIGATION  OF  A  CHEMICAL 
LASER  CAVITY  FLOWFIELD 


I.  Introduction 


1.1  Background 

Chemical  lasers  and  their  direct  application  to  the  Strategic  Defense  Initia¬ 
tive  (SDI  )  have  caused  a  large  increase  in  chemical  laser  research.  In  fact,  current 
concepts  for  the  Directed  Energy  Weapons  (DEW  )  program,  a  subprogram  of  SDI, 
have  stimulated  interest  in  four  promising  areas.  The  areas  are:  space-based  lasers 
(SBL  ),  ground-based  lasers  (GBL  ),  space-based  partical  beams  (SBPB  ),  and  nu¬ 
clear  directed  energy  weapons  (NDEW  ).  Decisions  on  whether  or  not  the  Chemical 
Laser  (CL  )  is  used  in  SDI  is  directly  related  to  the  performance  of  the  laser. 

Chemical  laser  performance,  or  more  specifically  the  beam  quality,  is  directly 
dependent  on  the  flow  structure  in  the  lasing  cavity.  Laser  performance  is  strongly 
connected  to  the  fluid  mechanics  processes  of  the  flow  crossed  by  the  beam.  Rapag- 
nani  and  Lankford  from  the  Air  Force  Weapons  Laboratory  (AFWL  )  state: 


Computer  modelling  of  the  flow  through  cavities  is  not  a  trival  task. 
Anaylsis  of  the  flows  are  complicated  by  the  existence  of  strong  cross¬ 
stream  pressure  gradients  and  subsonic  recirculation  zones  that  are  em¬ 
bedded  in  the  supersonic  cavity.  This  analysis  of  a  high  base  relief  cavity 
has  indicated  the  existence  of  large  property  variations  and  significant 
regions  of  subsonic  flow.  These  subsonic  recirculating  regions  are  essen¬ 
tially  dead  spaces  in  the  cavity  since  mass  is  neither  convected  in  or  out 
of  an  area  surrounded  by  a  closed  streamline.  These  regions  are  local  hot 
spots  since  the  temperature  is  approaching  its  stagnation  value.  Large 
pressure  variations  occurred  in  this  cavity.  (1:343-357) 
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Figure  1.1.  Schematic  of  a  Chemical  Laser  Nozzle  Assembly  and  Laser  Cavity 


Computer  analysis,  when  supported  by  experimental  results,  should  deepen 
the  understanding  of  many  complex  flows  associated  with  chemical  lasers.  From 
these  statements,  it  is  apparent  that  a  great  need  for  the  experimental  investigation 
of  a  chemical  laser  cavity  flowfield  exists. 


/ .  2  Scope 

Given  the  general  schematic  for  a  CL,  shown  in  Figure  1.1,  two  areas  of  concern 
must  be  investigated.  The  first,  and  primary,  area  of  concern  is  describing  the 
flowfield  in  the  cavity  and  its  limitations  under  simulated  laser  operating  conditions. 
Cold  flow  analysis  techniques  were  used  in  conjunction  with  theoretical  models  to 
describe  the  cavity  flowfield.  The  second  concern  is  the  design  and  configuration 
of  the  test  apparatus  required  to  support  the  test  section.  Pressure-vacuum  wind 
tunnel  technology  was  used  to  establish  and  maintain  laser  cavity-like  conditions. 
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1.3  Objectives 

The  investigation  has  two  major  objectives: 

1.  To  establish  a  facility  capable  of  investigating  flow  conditions  in  a  simulated 
laser  nozzle  assembly  and  lasing  cavity. 

2.  Based  on  time  history  of  pressures  in  the  cavity  and  schlieren  photographs  of 
the  cavity  flowfield,  determine  the  operating  characteristics  of  the  nozzle/cavity 
system. 

The  first  objective  was  accomplished  by  modeling  the  nozzle  and  cavity  design 
after  current  chemical  laser  nozzle  and  cavity  designs,  reconfiguring  a  current  blow¬ 
down  wind  tunnel  to  permit  simultaneous  flow  input  from  two  sources,  designing 
and  fabricating  a  mixer  and  flow  control  system  capable  of  providing  a  mixture  of 
air  and  helium  at  various  mass  flow  ratios,  and  designing  and  building  a  supersonic 
diffuser  capable  of  maintaining  the  test  section  cavity  at  the  low  pressures  required 
to  start  and  maintain  the  required  flow  conditions. 

The  last  objective  was  completed  by  a  combined  schlieren  and  static  pressure 
analysis  of  the  nozzle  and  cavity  flowfield.  One-dimensional  analysis  was  used  to 
describe  the  flow  conditions  in  the  nozzle  region.  A  two-dimensional  base-pressure 
model  was  used  to  analyze  the  flow  conditions  in  the  base  region  of  the  nozzle. 
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II.  Theory 


Analysis  of  the  cavity  flow  conditions  as  a  whole  is  nearly  impossible  with  the 
complicated  flow  patterns  and  their  interaction  with  each  other.  For  simplicity,  the 
flowfield  is  broken  up  into  four  logical  regions,  the  nozzle,  the  base,  the  cavity,  and 
the  variable  supersonic  diffuser  region.  The  conditions  in  the  first  region  were  analy¬ 
sed  assuming  one-dimensional  steady,  isentropic  conditions.  Schlieren  photographs 
of  the  nozzle  region  showed  no  shock  patterns  or  boundary  layer  separation  that 
would  preclude  the  use  of  the  above  mentioned  assumptions.  In  the  base  region, 
however,  the  flow  patterns  start  to  become  very  complicated.  The  underexpanded 
flow  stream  exiting  the  nozzle  goes  through  a  Prandtl-Meyer  expansion,  and  the 
flow  is  turned  and  accelerated.  Jet  boundary  streamlines  are  formed  indicating  the 
boundaries  of  the  accelerated  flow  stream.  A  multinozzle  arrangement  with  a  wide 
base  region  between  nozzles  creates  an  interaction  between  the  flow  streams  so  that 
when  the  accelerated  flows  meet  wake  regions  are  created  along  with  oblique  shocks 
and  recompression  regions.  To  handle  the  analysis  of  this  complicated  flow  stream, 
the  Korst  Generalized  Two-Dimensional  B else- Pressure  Model  was  selected  (3:593- 
595).  The  cavity  region  flowfield  is  then  a  series  of  oblique  shocks  with  wake  and 
recompression  regions.  Static  pressure  measurements  and  schlieren  photographs  were 
used  to  describe  the  conditions  in  the  cavity  flowfield.  Lastly,  it  was  necessary  to 
understand  the  flow  mechanics  in  a  supersonic  diffuser  so  that  one  could  be  designed 
to  help  maintain  the  test  conditions  in  the  cavity  for  a  longer  period.  Standard 
two-dimensional  diffuser  design  procedures  described  by  Hermann  (4:136-169)were 
used  to  design  the  diffuser. 
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Figure  2.1.  Flow  Model  for  Two-Dimensional  Base-Pressure  Problem  (3:594) 


2. 1  Nozzle  Theory 

One-dimensional  converging-diverging  (CD  )  nozzle  theory  described  in  Shapiro 
(5:76-105)  is  used  to  describe  the  flow  conditions  in  the  nozzle.  More  specifically, 
isentropic  equations  are  used  to  calculate  the  conditions  inside  the  CD  nozzle. 


2.2  Cavity  Base  Pressure  Flowfield  Analysis 

Analysis  of  the  cavity  base  pressures  was  handled  by  application  of  the  gen¬ 
eralized  restricted  theory  to  a  closed  wake  region(3:593-599).  Figure  2.1  shows  the 
flow  model.  Using  the  generalized  back-step  method  (3:597-599)  for  closed  wakes 
and  the  flow  configuration  shown  in  Figure  2.2,  an  iterative  technique  is  used  to 
calculate  flow  conditions  near  the  base,  and  the  theoretical  and  experimental  results 
are  compared.  The  iterative  procedure  is  detailed  below: 

1.  Select  A/la  and  0\a.  Where  9\a  is  the  streamline  angle  for  M\a. 
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Figure  2.2.  Flow  Configuration  for  a  Generalized  Back  Step  Closed  Wake  Region 
(3:599) 

2.  Calculate  i/i  for  M\a  Where  V\  is  defined  as  the  Prandtl-Meyer  function  for 

Ml0. 

3.  Calculate  ured.  Knowing  and  i/j,  i/rej  can  then  be  calculated  using: 

Vrti  =  V\  -  0U  (2.1) 

4.  Calculate  Mred,'  Using  an  iterative  technique  Mred,  can  be  calculated  with  the 
Prandtl-Meyer  equation  (6)  shown  below: 

v  =  6-Sarctan  — — -^]'5  -  arctan  ( M 2  —  l)'5  (2.2) 

0 


where 


5.  Guess  an  M2a.  Once  Mredt  has  been  calculated,  the  problem  is  handled  just  like 
a  basic  back  step  problem.  M2a  is  selected  close  to  Mredt  and  the  calculations 
are  continued. 


6.  Calculate  C2a  using: 


C2a  =  [1  + 


2 


where  Crocco’s  number  is 


Umax 


(2.3) 


(2.4) 


7.  Calculate  < p,.  Using: 


1  -  CIS 


rvR  < p2dr) 

-00  1  —  C,ap2 


where  p  is  called  dimensionless  velocity 


(2.5) 


<P  =  u/m«  (2.6) 

An  iteration  scheme  is  applied  solving  for  r),  where: 

V;  =  ^(l  +  er/77)  (2.7) 

and  t)r  =  2.  The  iteration  technique  described  above  was  performed  and  a 
polynomial  solution  for  as  a  function  of  C2a  was  obtained.  The  polynomial 
equation  is  listed  below: 

p,  =  3.15648  Cfa  -  9.22972  C|a  +  11.047  C24a 
-6.7842  Cfa  +  2.31116  C\a  -  0.382329  C2a  +  0.642023  (2.8) 
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8.  Calculate  C3y  Combine  Equations  2.4  and  2.6  to  form 


II 

J 

D 

(2.9) 

Since  C3a  =  C2a  then 

C2}  =  V3jC3a 

(2.10) 

implying  that 

^3]  =  f]C2a 

(2.11) 

Calculate  M3j  using: 

=  [ (CJ  -  1)1^V 

(2.12) 

10.  Calculate  From  the  escape  criterion  (3:596),  it  is  known  that  fluid  that  has 
^  <  1  will  not  be  able  to  penetrate  into  region  4  and  fluids  where  ^  >  1 
will  be  able  to  pentrate  region  4.  With  that  in  mind  and  the  fact  that  a  closed 
wake  specifies  that  SS2j-  —  ^  =  1 

r  P03d  Pi 

m  =(1  +  2T^Mi)^T  (2-13) 

P3  * 


11.  Calculate  02a.  From  jjj-,  M2a,  and  the  oblique  shock  relationships  (03a  —  04o) 
can  be  solved  for  using  (7:11): 


#3<i  —  #4«  =  arctan 


tt  Z  ~  1  .2,27Mia-(7-  l)-(7  +  l)^,.s 

lV7*?-f  +  lM  (7  +  lK  +  (7-l)  J 


(2.14) 


where 


Since  02a  =  03a  for  a  straight  jet  boundary  then 


Ola  —  0ia 


=  0ja  —  04a 
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12.  Calculate  vx  using: 

v\  —  —  <$12  (2.15) 

13.  Calculate  A iTede-  Given  Af2a,  <512,  and  vx,  A/redccan  be  calculated  using  Equa¬ 
tion  2.2  in  a  bi-section  iterative  technique.  The  resulting  Mredc  >s  then  com¬ 
pared  with  the  earlier  determined  Mredt ,  and  if  the  two  do  not  match,  a  new 
A/2a  is  selected  and  the  process  is  repeated.  The  process  is  continued  until  the 
calculated  Mredc  matches  Mred,. 

14.  Calculate  Knowing  A/red,  and  A/2a  and  the  isentropic  pressure  relationship 
listed  below  the  ratio  of  can  be  calculated  by: 

Pa  J 

£  =  (1  +  ~-V)^  (2.16) 

Po  2 

Since  the  flow  is  assumed  to  be  isentropic,  the  stagnation  pressure  p0  is  con¬ 
stant.  Therefore,  22  can  be  calculated  by: 

7  pi  J 

P2  _  Pi  P01 
Pi  P01  Pi 

15.  Compare  theoretical  and  experimental  The  calculated  &  can  now  be  com¬ 
pared  to  the  experimental  results.  At  this  particular  point,  it  should  be  noted 
that  theoretical  results  are  very  sensitive  to  the  original  9Xa  selected.  If  9Xa  is 
not  known  accurately  then  it  might  be  necessary  to  adjust  9Xa. 

2.3  Supersonic  Diffuser  Theory 

It  is  necessary  to  understand  the  shock  patterns  in  a  supersonic  diffuser  to 
determine  the  best  wedge  angle  for  the  convergent  section  of  the  diffuser.  The 
correct  shock  pattern  will  provide  the  best  (higher)  pressure  recovery.  To 
determine  the  correct  shock  pattern  some  of  the  diffuser  geometry  must  be 
known  and  an  optimization  must  be  performed  to  calculate  the  wedge  angle 


(2.17) 
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that  provides  the  best  recovery  pressure  for  the  diffuser  geometry  selected. 
Diffuser  theory  was  taken  from  Hermann  (4:125-136).  Appendix  B  lists  the 
software  program  and  some  of  the  calculations  necessary  in  the  selection  of  a 
diffuser  blade  wedge  angle. 
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III.  Experimental  Apparatus 


The  experimental  apparatus  is  composed  of  three  major  components:  the  flow 
system,  test  section  and  the  data  acquisition/reduction  system.  A  diagram  of 
the  system  is  shown  in  Figure  3.1.  The  flow  system  included  the  AFIT  com¬ 
pressed  air  system,  a  horizontal  stilling  chamber,  mixer,  a  variable  supersonic 
diffuser,  and  the  vacuum  system.  The  new  test  section  was  built  with  one 
set  of  hypersonic  nozzles  and  the  side  walls  were  made  of  optical  plexiglass  to 
allow  the  use  of  the  schlieren  flow  visualization  system.  Instrumentation  used 
to  monitor  the  flow  system  was  a  Qua  Tech  multichannel  simultaneous  sam¬ 
pling  system  installed  in  the  Zenith  Z-248  computer  and  ENDEVCO  pressure 
transducers  .  It  should  be  mentioned  that  a  flow  system  using  helium  and  air 
was  initially  considered,  but  it  was  finally  decided  to  use  solely  air  during  this 
research. 

3.0. 1  Flow  System  Considerable  design  and  fabrication  work  was  required  to 
set  up  the  new  flow  system  for  the  blowdown  wind  tunnel.  The  newly  designed 
mixer,  designed  and  fabricated  earlier  when  combined  air  and  helium  flows 
were  being  considered,  was  left  in  place  in  the  horizontal  stilling  chamber.  A 
supersonic  diffuser  was  designed  and  fabricated  to  attempt  to  maintain  the 
test  section  cavity  pressures  for  an  extended  time,  and  finally,  an  extension  to 
the  existing  vacuum  system  was  built  to  finish  the  flow  system. 

3.0.1. 1  Compressed  Air  System  The  compressed  air  system  consists  of  two 
ATLAS  COPCO  Air  Compressors,  Model  GAU-807,  two  DELTECH  Filters, 
Model  819,  two  Pioneer  Refrigerant  Air  Dryers,  Model  R550A,  and  two  Arrow 
Compressed  Air  Dryers.  Figure  3.2  shows  the  compressed  air  system  con¬ 
figuration.  Each  of  the  air  cooled  compressors  is  capable  of  providing  520 
standard  cubic  feet/minute  (SCFM  )  at  102  pounds/square  inch  gauge  (psig). 
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Figure  3.1.  Equipment  Apparatus  Configuration 


Figure  3.2.  Building  C40  Laboratory  Compressed  Air  System 


The  DELTECH  filters  are  used  to  eliminate  any  pollutants,  such  as  oil,  from 
the  air  flow.  The  Pioneer  refrigerant  air  dryers  eliminate  moisture  from  the 
air  flow  and  support  a  mass  flow  rate  up  to  500  SCFM.  The  calculated  pres¬ 
sure  loss  across  the  refrigerant  air  dryers  is  approximately  3.6  psig.  A  second 
drying  stage  is  available  through  the  Arrow  Chemical  Air  Dryers.  Normally, 
the  chemical  air  dryers  are  by-passed  due  to  their  large  pressure  loss,  approx¬ 
imately  10  psig.  Once  the  air  flow  is  through  the  final  stages  of  drying,  it  is 
piped  to  the  laboratories.  In  the  laboratory  the  air  goes  through  a  gate  valve 
and  into  a  cyclone  separator.  The  cyclone  separator  is  the  final  station  used  to 
eliminate  particulate  matter  from  the  air  flow.  Next,  the  air  is  routed  through 
an  electrically  controlled  3  inch  ball  valve.  The  Starlite  model  666-F  Ball  Valve 
provides  quick  opening  and  the  3  inch  ID  provides  the  maximum  air  flow  to 
fill  the  horizontal  stilling  chamber  and  start  the  nozzles. 
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Figure  3.3.  Mixer  End  Plate  View 

3. 0.1.2  Horizontal  Stilling  Chamber  The  six  foot  horizontal  stilling  chamber 
is  composed  of  two  parts.  The  forward  section  houses  a  steel  honeycomb  mesh 
to  help  straighten  and  distribute  the  air  evenly  prior  to  it  entering  the  mixer. 
The  forward  section,  shown  in  Figure  3.1  also  contains  the  mixer.  The  mixer 
provides  no  function  in  this  investigation.  The  second  section  of  the  horizontal 
stilling  chamber  acts  strictly  as  an  air  stilling  chamber  and  contains  no  other 
devices  other  than  pressure  transducer  taps. 

3. 0.1. 3  Mixer  The  mixer  design  was  based  on  the  Vassilatos  parallel  multi¬ 
tube  ejector  design  described  by  Zakanycz  (8).  The  dimensions  of  the  discharge 
end  of  the  Vassilatos  multitube  ejector  were  enlarged  so  that  a  large  scale  ver¬ 
sion  of  the  multitube  ejector  mixer  could  be  placed  inside  the  horizontal  stilling 
chamber  (see  Figure  3.3).  The  mixer  was  made  out  of  stainless  steel  tubing, 
plywood,  plexiglass,  8  inch  PVC,  and  aluminum  plate.  A  side  view  of  the  mixer 
design  is  shown  in  Figure  3.4.  The  mixer  slipped  into  the  forward  section  of 
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Air/Heliui  Mixer 


Figure  3.4.  Mixer  Side  View 

the  horizontal  stilling  chamber.  The  flange  end  of  the  mixer  was  sandwiched 
between  the  flange  ends  of  the  forward  and  rearward  portions  of  the  horizontal 
stilling  chamber.  Air  flow  entered  the  forward  portion  of  the  horizontal  stilling 
chamber  from  the  end  as  indicated  in  Figure  3.1.  The  air  flow  first  encountered 
a  steel  honeycomb  mesh  and  then  entered  the  mixer  through  the  forward  plate. 
Entering  the  mixer,  the  air  was  forced  to  flow  around  a  chamber  full  of  glass 
beads  (marbles).  The  glass  beads  caused  an  even  distribution  of  air  to  flow 
into  the  air  tubes.  A  fine  wire  mess  was  placed  between  the  glass  beads  and 
the  entrance  to  the  tubes  to  filter  out  any  debris.  Once  in  the  tubes,  the  air 
was  directed  through  the  remaining  part  of  the  mixer  and  out  the  end  plate. 

The  mixer  was  designed  to  receive  flow  from  a  second  source  which  would 
enter  a  cavity  formed  by  the  fit  of  the  mixer  installed  in  the  horizontal  stilling 
chamber.  The  side  wall  of  the  mixer  would  deflect  the  second  flow  stream 
around  the  mixer  and  into  the  holes  drilled  in  the  side  of  the  PVC.  Once  inside 
this  cavity  in  the  mixer,  the  second  flow  was  forced  to  flow  around  the  air 
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Figure  3.5.  Test  Section,  Top  View 


tubes  traversing  the  cavity,  into  a  second  set  of  tubes,  and  then  directed  out 
the  end  plate.  The  alternating  of  the  primary /secondary  tubes  exiting  the  end 
plate  would  induce  a  thorough  mixing  in  the  rearward  portion  of  the  horizontal 
stilling  chamber.  As  indicated  above,  a  second  flow  source  was  not  used  in  this 
investigation. 

3.0. 1.4  Variable  Supersonic  Diffuser  A  supersonic  diffuser  was  built  to  main¬ 
tain  the  low  pressures  in  the  cavity  portion  of  the  test  section.  The  diffuser  was 
made  out  of  aluminum  and  plexiglass  and  actuated  by  a  24  v  dc  Delco  motor. 
Figure  3.5  displays  the  mechanisms  required  to  vary  the  diffuser  dimensions. 
The  blades  of  the  diffuser  open  perpendicular  to  the  direction  of  expansion 
from  the  nozzles.  This  unique  design  was  adopted  to  minimize  the  length  of 
the  diffuser  section.  The  dc  motor,  see  Figure  3.5  was  electrically  controlled 
from  the  operator  control  panel.  Power  for  the  motor  was  provided  by  the 
laboratory  dc  power  supply,  a  Rapid  Electric  Co.  Model  S-528.  The  motor, 
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shown  in  Figure  3.6,  mounted  on  the  outside  of  the  diffuser  enclosure,  drives  a 
shaft  that  penetrates  the  diffuser  side  wall  through  a  pressure  tight  seal.  The 
shaft  drives  a  gearing  mechanism  that  turns  a  threaded  shaft.  The  threaded 
shaft  moves  back  and  forth  either  pushing  or  pulling  the  diffuser  slide  bracket. 
Adjustment  contacts  positioned  on  the  thread  shaft  deactivate  the  motor  when 
the  contact  bars  close  the  micro-switches  mounted  on  the  adjustment  bracket 
behind  the  gearing  mechanisms. 

3. 0.1. 5  Vacuum  System  To  establish  the  desired  test  section  operating  con¬ 
ditions  it  was  necessary  to  have  high  pressure  on  one  end  of  the  wind  tunnel 
and  a  vacuum  on  the  other  end.  The  vacuum  system  consists  of  16  large  tanks 
suspended  from  the  laboratory  and  three  vacuum  pumps.  Two  of  the  three 
vacuum  pumps  are  Leiman,  no  model  number,  radiator  cooled  pumps,  and  the 
third  pump  is  a  Stokes,  Model  212-H.  The  Stokes  pump  has  the  capacity  to 
pump  140  SCFM.  The  16  vacuum  tanks,  along  with  the  piping,  provide  ap¬ 
proximately  650  ft 3  of  vacuum  space.  A  network  of  8  inch  PVC  pipe  connected 
the  test  section  to  the  vacuum  system.  The  interface  between  the  PVC  pipe 
and  the  test  section  was  a  Consolidated  Vacuum  Corporation,  Model  VCS-41, 
high  vacuum  valve.  The  high  vacuum  valve  allowed  the  operator  to  isolate  the 
wind  tunnel  from  the  vacuum  system.  A  Meriam  U-tube  mercury  manometer, 
Model  20DA40,  was  incorporated  to  provide  the  operator  with  the  vacuum 
pressure  readings  while  the  tunnel  is  being  evacuated. 

3.0.2  Test  Section  The  test  section  was  designed  to  simulate  current  chemical 
laser  primary  nozzle  and  cavity  designs.  The  test  section  was  constructed  out 
of  aluminum  and  the  windows  were  cut  from  1  inch  thick  optical  plexiglass. 
Figure  3.7  shows  the  nozzle  assembly /cavity /diffuser  arrangement  of  the  test 
section. 


3-7 


3-8 


e  3.6.  Variable  Supersonic  Diffuser,  Top  View 


Figure  3.7.  Test  Section,  Side  View 


3.0.2. 1  Nozzles  The  nozzle/cavity  design  selected  and  its  general  operat¬ 
ing  conditions  were  established  from  information  provided  by  the  Air  Force 
Weapons  Laboratory  (  AFWL  )  (9),  Kirtland,  NM.  Figure  3.8  shows  the  con¬ 
tour  of  the  nozzle  and  Table  3.1  (10)  lists  the  dimensions  of  the  nozzle  contour. 

The  nozzle  design  parameters  are  listed  in  Table  3.2.  Figure  3.9  displays  the 
three  primary  wedges  as  they  are  mounted  in  the  test  section.  The  current 
design  allows  enough  material  to  machine  the  base  region  to  accommodate 
modifications  for  future  testing.  Three  pressure  ports,  shown  in  Figure  3.9 
were  drilled  into  each  of  the  primary  wedges  to  allow  static  pressure  measure¬ 
ments  in  the  region  of  the  base  and  the  nozzle  exit  plane. 


3.0.3  Data  Acquisition  and  Reduction  System  The  data  acquisition  system  is 
a  group  of  interconnected  electronic  instruments  used  to  measure  and  evaluate 
pressure  variables  produced  in  the  flow  stream.  The  data  reduction  system 
is  a  software  controlled  digital  computer  that  processes  and  analyzes  the  col- 
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Figure  3.8.  Nozzle  Contour 


Table  3.1.  Nozzle  Contour  Position 


Table  3;2.  Nozzle  Design  Parameters 


Variable 

Quantity 

7 

1.51 

Po 

100.0  psia 

To 

73  F 

Me 

6.00 

Pe 

0.104  psia 

Tt 

7.2  F 

AJA- 

28.8 

Figure  3.9.  Primary  Wedge  Configuration 
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Figure  3.10.  Data  Acquisition  System 


lected  data.  The  Data  Acquisition  and  Reduction  System  DARS  is  made 
up  of  several  basic  components.  A  block  diagram  of  this  system  is  shown 
in  Figure  3.10.  The  Zenith  Computer  System,  the  heart  of  the  DARS,  con¬ 
trols  the  data  acquisition  and  reduction  process  through  a  combined  Qua  Tech 
and  operator  written  Quick  Basic  software  program.  The  analog  signal  from 
the  piezoresistive  pressure  transducer  is  passed  to  the  Qua  Tech  Multichannel 
Simultaneous  Sampling  System  installed  in  the  Zenith  Computer  Processing 
Unit  CPU  .  The  Qua  Tech  sampling  system  allows  the  DARS  to  monitor  six 
channels  simultaneously. 

3.0.S.1  ENDEVCO  Piezoresistive  Pressure  Transducers  Miniature  piezore¬ 
sistive  pressure  transducers  with  highly  sensitive  strain  gauges  mounted  in  a 
silicon  diaphragm  were  used  to  monitor  the  flow  stream.  ENDEVCO  piezore¬ 
sistive  pressure  transducers  placed  in  the  flow  stream  converted  pressure  data 
to  electrical  analog  signals.  Pressure  measurements  were  taken  at  twenty-four 
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Table  3.3.  Pressure  Transducer  Signal  Conditioning  Unit  Assignments 


SCU 

Model  Number 

Serial  Number 

n 

8510B-5 

PP81 

I 

8510B-5 

79HB 

1 

8506B-5 

92BF 

fX 

8510B-5 

78HB 

1 

8506B-5 

74BF 

EX 

8506B-5 

68BF 

9 

8510B-15 

9BYP 

I 

8506B-5 

79BF 

ex 

8510B-100 

35DF 

positions  throughout  the  flow  system  and  the  test  section.  Four  different  sizes 
of  pressure  transducers  were  used.  One  ENDEVCO  pressure  transducer,  Model 
8510B-100,  was  ^  laced  in  the  side  wall  of  the  horizontal  stilling  chamber,  see 
Figure  3.1 1, upstream  of  the  test  section,  to  monitor  stagnation  pressure.  Three 
ENDEVCO  pressure  transducers,  Model  8510B-5,  were  used  in  the  nozzle  sec¬ 
tion  to  monitor  static  pressure,  and  four  ENDEVCO  pressure  transducers, 
Model  8506B-5,  were  used  in  the  test  section  plexiglass  side  walls  to  also  mon¬ 
itor  static  pressure.  One  ENDEVCO  pressure  transducer,  Model  8510B-15, 
was  placed  in  the  side  wall  of  the  vertical  stilling  chamber  to  monitor  the  back 
pressure.  Pressure  transducer  specifications  can  be  found  in  Reference  12.  The 
locations  of  the  twenty-four  pressure  taps  are  shown  in  Figure  3.11.  Since  a  six 
channel  data  acquisition  system  was  used,  only  combinations  of  six  positions 
could  be  measured  during  any  one  test  run.  Table  3.3  lists  the  model  and 
serial  number  of  the  transducers  used  on  a  particular  signal  conditioning  unit 
(SCU  ). 

3. 0.3. 2  ENDEVCO  Signaling  Conditioning  Units/Power  Supplies  Two  EN¬ 
DEVCO  Model  4225  ac-operated  power  supplies  were  used  to  provide  t  18 
volts  to  the  six  ENDEVCO  Model  4423  signal  Conditioning  units.  The  signal 
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Figure  3.11.  Pressure  Transducer  Locations 


conditioning  units  provide  excitation  voltage  to  the  pressure  transducer  and 
amplify  output  voltage  from  the  pressure  transducers.  Various  fixed  gain  set¬ 
tings  are  front  panel  selectable  and  were  used  to  match  the  input  voltage  range 
of  the  Qua  Tech  multichannel  simultaneous  sampling  system. 

3. 0.3. 3  Qua  Tech  Multichannel  Simultaneous  Sampling  System  The  Qua  Tech 
multichannel  system  provides  six  channels  of  simultaneous  sampling.  The 
six  channel  system  contains  three  PXB-721  parallel  expansion  boards,  six 
ADM12-11  12  bit  analog-to-  digital  (A/D  )  converter  modules,  one  CTM-11 
counter/timer  module  ,  and  one  SS-INT  interface  board.  An  input  voltage 
range  of  t  2.5  volts  was  selected.  The  amplified  analog  signal  was  sampled  at 
the  operator  specified  sample  rate  and  the  sampled  data  was  converted  to  an 
digital  format  by  the  A/D  card.  The  digital  signal  was  routed  to  the  computer 
CPU  over  the  computer  system  internal  wiring  system. 

3. 0.3. 4  Zenith  Computer  System  The  two  main  purposes  of  the  computer 
system  are  data  acquisition  and  data  reduction.  The  Zenith  Computer  System 
consists  of  one  Zenith  CPU,  Model  Z248,  and  one  Zenith  EGA  Color  Monitor, 
Model  ZVM-1380.  The  CPU  houses  a  20  M  hard  drive,  a  360  k  floppy  disk 
drive,  a  1.2  M  floppy  disk  drive,  a  math  co-processor  and  640  k  of  RAM. 
Qua  Tech  “LABSIM”  software  is  combined  with  an  operator  written  Quick 
Basic  software  program  to  control  data  acquisition  and  data  reduction.  The 
combined  software  package  and  its  use  is  shown  and  described  in  Appendix  A. 

3.0.4  Schlieren  Optical  System  Schlieren  photographic  equipment  was  used 
to  capture  the  flow  characteristics  in  the  test  section  cavity.  The  equipment 
configuration  is  shown  in  Figure  3.12  and  the  equipment  description  is  given 
in  Table  3.4. 
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Figure  3.12.  Schlieren  Equipment  Configuration 


Table  3.4.  Schlieren  Equipment  Data 


Item _ 

Spark  Lamp 

7.5  inch  Diameter  Concave  Mirrors 
Spark  Lamp  Power  Supply 
Trigger  D.C.  Power  Supply 
Knife  Edge 
Polaroid  Camera 
Polaroid  Film  Holder 


Model  Number 
Cordin  Model  5401 

Cordin  Model  5205 
Hewlett  Packard  (HP  ) Model  6205B 

Polaroid  Model  545 
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IV.  Experimental  Procedures 

Procedures  established  to  ensure  accurate  data  collection  cover  calibration  of 
instrumentation,  data  acquisition  and  reduction  procedures,  and  Schlieren  pho¬ 
tograph  procedures. 

4-1  Pressure  Transducer  Calibration 

AMETEK,  Pnuematic  Dead  Weight  Testers,  Model  PK-2  (low  pressure,  300 
inches  of  water)  and  Model  HK-500  (high  pressure,  500  psig),  were  used  to 
calibrate  the  ENDEVCO  piezoresistive  pressure  transducers.  The  test  configu¬ 
ration  is  shown  in  Figure  4.1.  The  manufacturer  of  the  pnuematic  dead  weight 
testers  cautions  that  the  use  of  any  gas  with  the  dead  weight  testers  other  than 
clean  compressed  air  or  nitrogen  could  cause  discrepancies  in  the  calibration 
data  (13).  The  pressure  transducers  were  calibrated  with  the  ENDEVCO  SCU 
in  line.  The  ENDEVCO  SCU  has  four  gain  settings  5,  10,  20,  50  and  a  bridge 
balance  adjustment  screw.  Since  the  ENDEVCO  SCU  does  not  have  a  unity 
gain  setting,  it  will  always  multiply  the  pressure  transducer  output  voltage 
by  some  factor  depending  on  which  gain  setting  has  been  selected.  The  gain 
factors  were  not  exactly  as  indicated  on  the  front  panel;  therefore,  the  pressure 
transducers  were  calibrated  with  the  SCU  in  line  so  that  the  exact  gain  factor 
is  included  in  the  slope  of  the  calibration  curve.  The  calibration  curves  of  the 
pressure  transducers  used  in  this  research  are  presented  in  Appendix  B. 

4-2  Data  Collection 

Two  methods  of  data  collection  were  used:  data  acquisition  and  reduction  of 
pressure  measurements,  and  schleiren  photographs  of  the  flowfield. 

4-2.1  Data  Acquisition  and  Reduction  Procedures  Data  acquisition  and  re¬ 
duction  is  controlled  by  the  operator  software  described  in  Appendix  B.  The 
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Figure  4.1.  Pressure  Transducer  Calibration  Setup 


software  has  three  distinct  sections.  The  first  section  handles  data  acquisition. 
The  second  section  handles  data  reduction,  and  the  third  provides  a  visual  look 
at  the  converted  data  on  the  monitor  screen.  Each  section  walks  the  operator 
through  the  procedures  necessary  to  complete  the  operation.  See  Appendix  B 
for  more  details. 

4.2.2  Schlieren  Photograph  Procedures  Schlieren  photographs  were  used  to 
record  the  visual  representation  of  the  flow  conditions  in  the  test  section  cavity. 
Photographs  were  taken  at  specific  times  during  the  test  runs.  Pictures  were 
then  compared  to  the  flow  data  collected  by  the  DARS.  The  actual  procedures 
for  setting  up  and  operating  schleiren  photographic  equipment  are  described 
in  AGARDograph  No  23  (14:13-15). 
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V.  Results  and  Discussions 


The  flow  system  consisted  of  a  primary  source  of  air  fed  into  a  pressure- vacuum 
wind  tunnel  test  section  through  a  mixing  device.  Once  through  le  test 
section  the  air  flow  was  routed  through  a  supersonic  diffuser,  and  emptied  into 
a  vacuum  system.  System  evaluation  included  vacuum  performance,  diffuser 
performance,  and  test  section  (nozzle/base/cavity)  performance. 

5.1  Flow  System  Equipment  Performance 

The  complete  installation  with  controls  and  instrumentation  was  tested  to 
determine  its  suitability  to  simulate  the  necessary  laser  nozzle/cavity  operation. 
The  system  was  able  to  pull  a  vacuum  down  to  0.1  psia.  With  the  vacuum 
pumps  turned  off  and  the  vacuum  system  closed,  the  vacuum  was  maintained 
at  0.1  psia  for  approximately  1  hour. 

The  variable  supersonic  diffuser  needed  to  maintain  the  low  cavity  pressures 
was  designed  using  two-dimensional  theory.  Diffuser  theory  indicated  that 
the  optimum  starting  and  minimum  running  throat  area  ratios  for  a  diffuser 
with  an  entrance  Mach  number  of  5.0,  should  be  0.6  and  0.2,  respectively. 
Experimentally,  the  optimum  starting  and  minimum  running  throat  area  ratios 
were  never  established.  Data  runs  indicate  that  the  minimum  running  area 
ratio  was  much  higher  than  the  theoretical  value.  The  minimum  diffuser  throat 
area  ratio  that  allowed  the  full  starting  of  the  nozzles  and  held  a  steady  static 
pressure  at  transducer  positions  17  and  18  was  determined  to  be  0.625.  Diffuser 
area  ratios  above  0.625  also  allowed  the  starting  of  the  nozzles  but  failed  to 
maintain  the  lower  pressures  for  as  long  as  the  0.625  diffuser  throat  setting. 
Throat  areas  below  the  0.5625  setting  do  not  allow  the  nozzle  flow  to  fully 
expand  in  the  cavity  and  in  some  cases  cause  separation  in  the  nozzle  divergent 
section. 
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It  was  initially  intended  that  the  flow  testing  would  be  done  with  an  air/helium 
mixer,  but  it  was  later  decided  to  do  the  initial  testing  using  strictly  air.  Even 
though  the  mixer  was  not  used  for  its  intended  purposes,  it  was  still  placed 
in  the  flow  stream.  Pressure  across  the  mixer  was  monitored  and  no  pressure 
drop  was  observed.  Once  the  major  components  were  in  place,  an  electrically 
controlled  ball  valve  was  used  to  quickly  turn  the  air  flow  on  and  off.  Testing 
indicated  that  the  horizontal  stilling  chamber  reached  maximum  pressure  in 
approximately  2  seconds. 

5.2  Test  Section  Performance 

The  test  section  cavity  flowfield  conditions  were  investigated  using  two  tech¬ 
niques,  schlieren  photography  and  static  pressures  measurements.  Photographs, 
as  well  as  pressure  measurements,  were  taken  with  various  supersonic  diffuser 
throat  areas.  Examples  of  two  of  the  test  flow  runs  are  shown  in  Figures  5.1 
and  5.2.  Close  inspection  of  the  photographs  will  show  a  multitude  of  fluid 
dynamics  interaction  such  as  boundary  layers  by  the  side  walls  of  the  noz¬ 
zles,  wakes  and  shock  reflections  along  the  base  centerline,  and  jet  streamlines 
extending  from  the  exit  of  the  nozzle  to  the  base  centerline.  All  the  flow  condi¬ 
tions  mentioned,  and  even  some  that  were  not  mentioned,  have  an  effect  on  the 
accuracy  of  the  experimental  data.  For  this  reason,  the  analysis  was  broken 
up  into  three  areas,  the  nozzle,  the  base,  and  the  open  cavity  section.  Detailed 
investigation  of  the  flow  in  the  diffuser  was  not  accomplished. 

To  start  the  investigation  of  the  flowfield  in  the  cavity,  a  series  of  schlieren 
photographs  were  taken.  The  photographs  were  used  to  select  locations  for 
pressure  transducers  in  the  test  section  plexiglass  sidewall  based  on  observed 
wakes  and  shock  wave  positions.  Pressure  transducers  were  also  located  in 
the  horizontal  stilling  chamber  and  downstream  from  the  diffuser  to  provide 
both  static  and  stagnation  pressure  readings.  Nine  other  pressure  transducer 
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Figure  5.1.  Cavity  Flowfield,  Diffuser  Throat  Setting  of  0.6875  inches 
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Figure  5.3.  Cavity  Flowflield,  Diffuser  throat  setting  of  0.75  inches 


taps  were  built  into  the  nozzle  design,  three  taps  per  nozzle,  to  provide  addi¬ 
tional  static  pressure  readings.  Figure  3.11  provides  a  complete  picture  of  the 
exact  pressure  transducer  locations.  Figures  5.3-5. 6  show  the  flowfields  with 
the  supersonic  difusser  throat  set  at  0.75,  0.6875,  0.625,  and  0.5625  inches, 
respectively. 

5.2.1  Nozzle  Flowfield  Analysis  The  design  parameters  of  the  nozzle  are 
listed  in  Table  3.2.  From  Figures  5. 3-5. 6,  it  can  be  seen  that  the  exit  con¬ 
ditions  from  the  suppersonic  nozzle  seem  to  be  shock  free.  For  this  reason,  the 
analysis  of  the  nozzle  was  handled  using  general  one-dimensional,  steady,  isen- 
tropic  flow.  The  nozzle  was  designed  to  operate  with  a  7  of  1.51,  but  this  intial 
series  of  testing  used  strictly  air,  so  it  was  necessary  to  recalculate  the  exit 
Mach  number  (Af*  )of  the  nozzle.  Using  the  Ae/A*  of  28.8  and  a  7  of  1.4  the 
new  Me  was  calculated  to  be  5.18.  The  horizontal  stilling  chamber  sidewall 
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Figure  6.5.  Cavity  Flowflield,  Diffuser  throat  setting  of  0.625  inches 


Figure  5.6.  Cavity  Flowflield,  Diffuser  throat  setting  of  0.5625  inches 
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pressure,  taken  prior  to  the  test  section,  was  considered  to  be  a  stagnation 
pressure  due  to  the  low  velocities  in  the  horizontal  stilling  chamber.  Using 
the  static  pressure  taps  located  at  the  exit  of  the  nozzles,  see  Figures  5. 3-5. 6, 
and  the  stagnation  pressure,  the  exit  Mach  number  for  all  four  diffuser  throat 
settings  was  calculated  and  is  shown  in  Figures  5.7-5.10.  The  exit  Mach  num¬ 
bers  are  all  very  close  to  the  theoretical  value  of  5.18  which  indicates  that 
the  assumption  of  one-dimensional,  steady,  isentropic  flow  was  correct.  The 
Mach  number  falloff  shown  in  Figures  5.9-5.10  is  caused  by  the  increase  in 
static  pressure  at  the  nozzle  exit  plane,  which  is  brought  on  by  the  less  efficent 
diffuser  throat  area  setting. 

5.2.2  Base  Flow  Region  Analysis  Analysis  of  the  flowfield  after  the  nozzle 
exit  is  more  complicated  than  the  one-dimensional  analysis  used  in  the  nozzle 
section.  Figures  5. 3-5. 6  show  various  flow  conditions  depending  on  the  value 
of  the  diffuser  throat  area.  Close  to  the  exit,  though,  and  in  the  base  region 
it  is  possible  to  use  the  two-dimensional  Korst  theory  (3:593-599)  to  calculate 
base  pressures  and  Mach  numbers  further  down  stream  of  the  exit.  To  use 
the  Korst  theory,  a  M\a  and  0i„  must  be  selected.  The  experimental  Mach 
number  calculated  by  the  static  pressure  at  the  exit  of  the  nozzle  was  used  as 
M\a.  The  flow  stream  angle,  however,  is  not  quite  as  easy.  From  the  method  of 
characteristics,  it  can  be  shown  that  the  flow  stream  angle  varies  from  zero  to 
the  maximum  deflection  starting  at  the  centerline  and  moving  out  to  the  edge 
of  the  nozzle  (see  Figure  5.11).  Additionally,  the  Mach  number  varies  from 
a  minimum  to  a  maximum  starting  at  the  nozzle  centerline.  To  compare  the 
experimental  results  with  a  theoretical  value  it  was  first  necessary  to  determine 
the  theoretical  mass  averaged  flow  stream  angle.  Since  the  divergent  section 
of  the  nozzle  has  a  slope  of  approximately  twenty-six  degrees,  the  flow  stream 
angle  should  be  somewhere  between  0  and  26  degrees.  Since  the  Mach  number 
increases  towards  the  outside  edge  of  the  nozzle,  the  mass  flow  rate  would 
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Figure  5.7.  Nozzle  Exit  Mach  Number  Me  vs  Time,  Diffuser  Throat  0.5625  inches 
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Figure  5.9.  Nozzle  Exit  Mach  Number  Me  vs  Time,  Diffuser  Throat  0.6875  inches 
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Figure  5.10.  Nozzle  Exit  Mach  Number  Me  vs  Time,  Diffuser  Throat  0.75  inches 
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Figure  5.11.  Variance  of  Mach  Number  and  Flow  Stream  Angle  at  the  Exit  Plane 
of  a  Supersonic  Nozzle 


also  increase  causing  the  center  of  mass  to  shift  towards  the  outside  edge. 
Shifting  the  center  of  mass  to  the  outside  edge  would  also  shift  the  average 
flow  angle  towards  the  outside  edge.  With  that  in  mind,  the  Korst  method 
was  then  used  to  determine  the  flow  angle  that  fit  the  experimental  pressure 
The  average  of  the  experimentally  determined  flow  angles  was  calculated, 
18.4  degrees,  and  was  used  as  ^ia,Htor„ical-  Using  the  experimentally  calculated 
Me  and  &\ atH4er„itali  theoretical  values  for  the  pressures  and  the  Mach  numbers 
in  the  base  region  were  calculated.  The  difference  between  the  experimental 
pressure  ratio  and  the  theoretical  pressure  ratio  were  compared  to  indicate 
the  validity  of  using  the  two-dimensional  generalized  Korst  model.  Table  5.1 
summarizes  the  results  of  several  calculations.  Figure  5.12  shows  the  Korst 
theoretical  curve  and  a  data  point  calculated  using  the  experimental  data. 
Figure  5.13  demonstrates  the  effect  of  various  flow  stream  angles  as  well  as 
positions  one  of  the  experimentally  derived  points.  From  the  results  it  would 
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Table  5.1.  Base  Region  Results  Comparison 


Experimental 

Theoretical 

Difference 

$la 

El 

Pi 

e,a 

El 

pi 

(Percent) 

0.9012 

wm 

0.8893 

1.33 

0.8995 

EH 

0.8880 

1.29 

18.5 

0.9051 

m 

0.8924 

1.41 

18.2 

0.8622 

18.4 

0.8840 

2.50 

seem  that  the  Korst  method  is  an  accurate  model  of  the  base  region  flow 
conditions.  Further  investigation  into  the  selection  of  0ia,(Uor«„eoj  is  necessary 
to  confirm  the  accuracy  of  the  results. 

5.2.3  Cavity  Flowfield  Anaysis  Static  pressure  taps  and  schlieren  photographs 
were  used  to  analyze  the  flow  conditions  in  the  cavity.  Figures  5. 3-5. 6  show 
various  flow  conditions  at  different  diffuser  throat  settings.  The  first  three 
diffuser  throat  settings  produce  underexpanded  conditions  at  the  exit  of  the 
nozzles  and  the  flow  expands  outward  to  form  the  diamond  pattern  shown  in 
the  photographs.  If  the  flow  between  transducer  positions  17  and  14,  shown 
in  Figure  5.8,  is  considered  to  be  shock  free,  then  isentropic  conditions  can  be 
applied.  Using  the  stagnation  pressure  from  upstream,  the  Mach  number  at 
postion  14  can  be  calculated.  Figure  5.14  shows  that  the  Mach  number  at  pos¬ 
tion  14  is  initially  larger  than  the  Mach  number  at  position  17,  this  indicates 
that  normal  isentropic  expansion  applies  for  the  first  few  seconds  of  the  run 
and  then  cavity  conditions  begin  to  detoriate  producing  a  lower  Mach  num¬ 
ber.  Figures  5.15  and  5.16  show  that  static  pressure  increases  long  the  nozzle 
centerline,  with  respect  to  time.  Additionally,  static  pressure  at  position  11 
was  noted  to  jump  above  the  static  pressure  at  position  9,  possibly  indicating 
a  strong  shock  system  moving  past  position  11  and  back  towards  the  nozzle 
exit.  Further  investigation  is  necessary  to  understand  why  position  9  pressure 
readings  are  less  than  position  11.  Pressure  values  in  the  figures  where  the 
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Figure  5.13. 
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Effect  of  Approaching  Flow  Stream  Angle  in  Korst  Generalized  Back 
Step  Calculations,  Diffuser  Throat  Setting  0.625  inches 


phrase  smoothed  data  is  used  indicates  the  pressure  data  has  been  artifically 
smoothed  using  an  FFT  program.  Both  Figure  5.15  and  5.16  indicate  low, 
stable  static  pressure  conditions  close  to  the  exit  of  the  nozzle,  but  pressure 
becomes  more  unstable  further  out  in  the  cavity  section.  Figures  5.17-5.20 
show  the  base  centerline  static  pressures,  smoothed  data,  for  all  four  diffuser 
settings.  Trends  in  the  smoothed  data  indicated  the  effects  of  the  different 
diffuser  settings.  For  example,  Figure  5.18  shows  that  overall  static  pressure 
remains  lower  longer  with  the  0.625  setting  than  with  the  other  three  diffuser 
settings.  However,  analyzing  the  unsmoothed  data  points  also  provides  some 
insight  into  the  dynamics  of  the  flow  stream.  Figures  5.21-5.32  show  both 
the  actual  and  the  smoothed  datapoints  for  two  diffuser  settings,  0.625  and 
0.5625  and  clearly  indicate  the  fluxuating  conditions  in  the  cavity.  A  compari¬ 
son  of  the  figures  indicate  that  the  0.625  diffuser  setting  provides  stable  cavity 
conditions  for  the  longest  time  period. 
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Figure  5.14.  Nozzle  Centerline  Mach  Number  vs  Time,  Diffuser  Throat  Area  Ratio 
0.625 
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Figure  5.15.  Nozzle  Centerline  Static  Pressure  Measurements  vs  Time,  Diffuser 
Throat  Area  Ratio  of  0.625,  Smoothed  Data 
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Figure  5.16.  Nozzle  Centerline  Static  Pressure  Measurements  vs  Time,  Diffuser 
Throat  Area  Ratio  of  0.5625,  Smoothed  Data 
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Figure  5.17.  Base  Centerline  Static  Pressure  Measurements  vs  Time,  Diffuser 
Throat  Area  Ratio  of  0.5625,  Smoothed  Data 
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Figure  5.18.  Base  Centerline  Static  Pressure  Measurements  vs  Time,  Diffuser 
Throat  Area  Ratio  of  0.625,  Smoothed  Data 
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Figure  5.20.  Base  Centerline  Static  Pressure  Measurements  vs  Time,  Diffuser 
Throat  Area  Ratio  of  0.75,  Smoothed  Data 
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Figure  5.21.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.5625,  Trans¬ 
ducer  Position  #1 
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Figure  5.22.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.5625,  Trans¬ 
ducer  Position  #14 
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Figure  5.23.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.5625,  Trans¬ 
ducer  Position  #7 
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Figure  5.24.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.5625,  Trans¬ 
ducer  Position  #16 
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Figure  5.25.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.5625,  Trans¬ 
ducer  Position  #8 
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TIME,  SECONDS 


Figure  5.27.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.625,  Trans¬ 
ducer  Position  #1 
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Figure  5.28. 
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TIME,  SECONDS 


Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.625,  Trans¬ 
ducer  Position  #14 
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<  0.00  10.00  20.00 

TIME,  SECONDS 

Figure  5.29.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.625,  Trans¬ 
ducer  Position  #7 
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Figure  5.30.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.625,  Trans¬ 
ducer  Position  #16 
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Flow  conditions  in  the  cavity  regions  are  very  complex.  A  more  thorough 
investigation  into  the  individual  flow  charateristics  is  necessary  before  the  flow 
patterns  can  be  fully  understood.  Additional  information  such  as  stagnation 
pressure  in  the  cavity  and  additional  static  pressure  ports  might  improve  the 
mapping  of  the  flowfield.  Trends  in  the  base  centerline  pressures  can  also  help 
analyze  the  wake  closure  and  wake  opening.  Comparision  of  the  cavity  static 
pressures  with  the  vacuum  static  pressure,  transducer  position  5,  indicates 
that  the  base  centerline  wake  regions  start  to  open  up  when  the  centerline  to 
vacuum  static  pressure  ratio  p~v£*u  is  higher  than  0.20.  Further  investigation 
into  the  possible  application  of  the  Korst  escape  criterion  to  the  wake  might 
help  describe  the  flow.  Finally,  interferometer  photographs  might  improve 
visual  inspection  of  the  low  density  flowfield. 
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Figure  5.32.  Static  Pressure  vs  Time,  Diffuser  Throat  Area  Ratio  of  0.625,  Trans¬ 
ducer  Position  #6 


VI.  Conclusions 


An  experimental  investigation  of  the  flowfield  in  a  simulated  chemical  laser 
cavity  was  accomplished.  The  time  history  of  static  pressure  at  various  loca¬ 
tions  in  the  cavity  and  schlieren  photographs  of  the  entire  flowfield  formed  the 
basis  for  the  performance  evaluation. 

(a)  The  facility  established  to  investigate  the  flow  conditions  in  a  simulated 
laser  nozzle  assembly  and  lasing  cavity  is  suitable  for  this  type  of  cold 
flow  investigation. 

(b)  The  nozzles  performed  at  the  expected  design  conditions.  Design  condi¬ 
tions  are  established  in  the  nozzle  section  immediately  after  flow  initiation 
for  various  supersonic  difusser  throat  area  ratios.  Conditions  are  main¬ 
tained  for  various  lengths  of  time  and  are  dependent  on  the  build-up  of 
back  pressure  within  the  cavity.  Results  indicate  that  the  nozzle  and  base 
regions  followed  theory  very  closely. 

(c)  The  combined  Korst  general/basic  back  step  theory,  when  carefully  ap¬ 
plied,  accurately  predicts  experimental  base  pressure  ratios. 

(d)  Flow  conditions  in  the  cavity,  determined  by  static  pressure  taps  and 
schlieren  photographs,  are  directly  related  to  the  supersonic  diffuser  throat 
setting.  Settings  below  the  minimum  running  throat  area  ratio  caused  the 
flow  pattern  to  breakup,  opening  the  base  wake  regions,  and  unstarting 
the  nozzle  or  causing  the  separation  of  the  flow  stream  inside  the  nozzle 
exit. 

(e)  Limited  data  on  the  supersonic  diffuser  indicates  that  it  was  not  operating 
like  the  two-dimensional  model  used  to  design  it,  but  flow  conditions  in  the 
cavity  were  being  maintained  at  operation  conditions  longer  than  would 
be  possible  without  the  diffuser.  Three-dimensional  flow  patterns  are 
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suspected  to  be  the  main  reason  why  the  diffuser  did  not  operate  closer 
to  the  two-dimensional  model. 


6-2 


VII.  Recommendations 


The  study  of  flow  conditions  in  the  cavity  of  a  chemical  laser  is  extremely 
important  to  the  furthering  of  chemical  laser  research.  Proper  cavity  conditions 
are  a  must  to  insure  safe  and  efficient  lasing.  Cold  flow  analysis  of  the  cavity 
flowfield  provides  a  foundation  for  further  hot  flow  studies  and  will  provide  a 
clearer  understanding  of  complicated  flowfields.  The  following  subjects  require 
further  investigation: 

(a)  Install  total  pressure  probe  in  the  cavity  section  for  measurement  of  total 
pressure  and  correlation  with  static  pressure  values. 

(b)  Install  pressure  devices  in  the  diffuser  to  describe  and  analyize  the  flowfield 
conditions. 

(c)  Install  additional  static  pressure  ports  in  the  test  section  sidewall  to  more 
completely  map  out  the  flowfield. 

(d)  Install  and  test  an  air/helium  flow  control  system  that  actively  monitors 
and  regulates  the  air/heliurn  mixture  through  the  flow  system. 

(e)  Investigate  other  nozzle  configurations  that  allow  base  and  hypersonic 
wedge  flow. 
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Appendix  A.  Computer  Programs 


A. I  Diffuser  Program 

To  determine  the  diffuser  blade  geometry,  it  was  first  necessary  to  know  the 
entrance  Mach  number  to  the  diffuser,  the  wedge  angle  (u>  ),  the  specific  heat 
ratio  (7  ),  the  number  of  shocks  to  be  considered,  and  the  diffuser  throat 
area  ( Apt  )•  A  program  was  written  to  calculate  the  conditions  after  several 
oblique  shocks.  Using  the  diffuser  program,  the  optimum  recovery  pressure 
was  determined  for  various  bl?Je  geometries.  The  blade  geometries  were  then 
compared  to  the  acceptable  diffuser  overall  dimensions,  and  the  best  blade 
design  was  selected.  The  following  is  a  copy  of  the  Fortran  program  used  to 
calculate  some  of  the  diffuser  shock  conditions. 


C234567 

C 

C 

PROGRAM  DIFFUSER 
IMPLICIT  REAL+8  (A-Z) 

IITEGER  FLAG.I.I.IEG 

DIMEISIQI  A(4) ,B(4) ,C(4) ,FLAG(10) ,M(10) ,PS(10) ,PT(10) ,WSA(10) , 
1  XR(3) ,XI(3) 

C0MM0I  K. RAD, DEC 

SDEBUG 

C 

C 

C  *  INTERACTIVE  IIPUT  OF  KIOWIS.  * 

c 

WRITE(* , 10) 

10  FORXAKSX,  'EITER  MACH  HUMBER  1:  *) 

READ(*,*)MX 

WRITE(*,20) 

20  FORMAT (6X, 'EITER  VEDGE  AHGLE:  ') 

Rr  >,*)VAD 
MKaTE(*,30) 

30  F0RMAT(SX, 'EITER  THROAT  WIDTH  (IICHES):  ') 

READ(*,*)TW 
WRITE(* ,40) 
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40  F0RMAT(5X,'EITER  HUMBER  OF  OBLIQUE  SHOCKS  TO  BE  CONSIDERED:  ') 
READ(*,*)I 
WRITE(*,4S) 

45  FORMAT(SX,'EITER  GAMMA:  ') 

READ(* , *)K 
C 
C 

c  ****************************************************** 

C  *  SET  COISTAITS .  * 

C  ****************************************************** 

c 

DEG*5 . 729577951301 
RAD= 1 . 745329252D-02 
A1=9.27D0 
AT=TW*9.27D0 
C 
C 

c  ****************************************************** 

C  *  I1ITIALIZE  VALUES.  * 

C  ****************************************************** 

c 

WAR=WAD*RAD 

M(1)=MX 

PTT=1D0 

C 

C  ****************************************************** 

C  •  DETERMIIE  THE  ROOTS  OF  THE  POLYIOMIAL  USIIG  THE  * 

C  *  BAIRSTOW’S  METHOD.  THE  SMALLEST  OF  THE  THREE  * 

C  *  ROOTS  CORRESPOHDS  TO  A  DECREASE  II  EITROPY  * 

C  *  AID  SHOULD  BE  DISREGARDED.  THE  REMAIIIIG  • 

C  *  TWO  ROOTS  WILL  BE  THE  STROIG  AID  THE  WEAK  * 

C  *  SHOCK  WAVE  AIGLES .  * 

C  ****************************************************** 

c 

DO  50  1*1,1 

CALL  SAMAX(H(I) ,SAMR) 

CALL  WAMAX(SAMR,M(I),WAMR) 

IF(WAR.LE.WAMR)THEI 
CALL  PCIAR(WAR , M ( I ) , SS AR , WSAR) 

CALL  OBSHOC(WSAR,M(I),MY,PSYX,PTYX) 

M(I+1)*MT 

PS(I)*PSYX 

PT(I)*PTTX 

WSAD*WSAR*DEG 

WSA(I)*WSAD 

FLAG(I)*1 

EIDIF 

IF (WAR . GT . WAMR)THEI 
WRITE(* ,60)M(I) 

60  F0RMAT(5X, 'THE  WEDGE  AIGLE  EXCEDES  THE  MAXIMUM  WEDGE  AIGLE ’ , 

1  6X, *  FOR  MAC.  7:  ' .F5.2) 
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GOTO  1000 
EIDIF 

50  COITIIUE 

CALL  MFLUX(MCl).MFl) 

CALL  MFLUX(M(I) ,MFI) 

DO  70  1=1,1 

PTT=PTT*PT(I) 

70  COITIIUE 
AR=AT/A1 

AR1=MF1/(MFI*PTT) 

C 

C 

C  ****************************************************** 

C  *  PRUT  RESULTS.  * 

C  •  ****************************************************** 

c 

WRITE (* ,150)AR,AR1,M(I+1),K 
150  FORMAT (1 IX, 'TRUE  DIFFUSER  AREA  RATIO:  '.F6.3, 

1  /,5X, ’CALCULATED  DIFFUSER  AREA  RATIO:  ’,F6.3, 

2  /,5X, ’MACH  IUMBER  AFTER  LAST  OBLIQUE  SHOCK:  ’, 

3  F6 . 3 , / , 5X , ’ RATIO  OF  SPECIFIC  HEATS  (GAMMA):  *,F8.3) 

WRITE(* ,200) 

200  FORMAT (/, 5X , ’REGIOI’ ,5X, ’ MACH# »,9X,’ SHOCK  AIGLE’,7X,’  PR  STATIC’ 
1  91, ’PR  TOTAL’,/) 

DO  220  1=1,1 

IF(FLAG(I) .EQ.O)GOTO  1000 

WRITE(« ,210)1 ,M(I) , WSA(I) , PS(I) ,PT(I) 

210  F0RMAT(6X,I2,4X>F12.9.4X,F12.9,4X,F13.9,4X,F13.9) 

220  COITIIUE 
1000  STOP 
EID 
C 

SUBROUTIIE  MFLUX(M.MF) 

IMPLICIT  REAL*8  (A-Z) 

IITEGER  I.I.IEG 
COMMOI  K. RAD, DEG 
MF1=2D0/(K+1D0) 

MF2=((K-1DO)/2DO)*M**2DO 

MF3*-(K+1DO)/((K-1DO)*2DO) 

MF=M*(MF1*(1D0+MF2) )**MF3 
RSTURI 
EID 
C 
C 

c 
c 
o 
c 
c 
c 
c 


*  SUBROUTIIE  PCHAR 

*  THIS  SUBROUTIIE  CALCULATES  THE  CHARACTERISTICS 

*  OF  THE  POLYIOMIAL  EQUATIOI  AID  TBEI  CALLS  POLY 

*  TO  SOLVE  THE  EQUATIOI. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUT I IE  PCHAR(HAR, M , SSA , MSA) 

IMPLICIT  REAL*8  (A-Z) 

IITEGER  I.I.IEG 

DIMEISIOI  A(4) ,B(4) ,C(4) ,XR(3) ,XI(3) 

COMMOI  X, RAO, DEG 
1=3 

A(0)=-1D0*(((DC0S(WAR))**2D0)/M**4D0) 

A1A=( (2DO*M**2DO)+1DO)/M**4DO 
AlB=((K+l)**2DO)/4DO  +  (K-1)/M**2D0 
A(1)=A1A+A1B*(DSII(WAR))**2D0 

A(2)=-1DO*((M**2DO+2DO)/M**2DO)-(K*(DSII(WAR))**2DO) 

A(3)=1D0 

CALL  POLY (I , A , XR , XI , HEG) 

1=3 

DO  10  1=1,1 

IF(XR(I).LT.ODO)GOTO  10 
XR(I)=DSQRT(XR(I)) 

XR(I)=DASII(XR(I) ) 

XR(I)=XR(I)*DEG 
10  COITIIUE 

SSA=DMAX1(XR(1) ,XR(2) ,XR(3>) 

DO  20  1=1,1 

HI1=XR(I)/SSA 
IF(HI1 . EQ. 1)THEI 
XR(I)=0D0 
EIDIF 
20  COITIIUE 

VSA=DHAX1(XR(1) ,XR(2) ,XR(3)) 

SSA=SSA*RAD 

VSA=VSA*RAD 

RETURI 

EID 


*  SUBROUTIIE  POLY 

« 

*  THIS  SUBROUTIIE  SOLVES  THE  THIRD  ORDER  POLYIOMIAL 

*  FOR  THE  THREE  ROOTS.  TWO  OF  THE  ROO. '  ARE  THE 

*  STROIO  AID  WEAK  SHOCK  AIGLES.  THE  THIRD  ROOT 

«  THE  SMALLEST  VALUE  CORRESPOIDS  TO  A  DECREASE  II 

*  EITROPT  AID  SHOULD  BE  DISREGARDED. 


SUBROUTIIE  POLY(I,A,XR,XI,IEG) 
IMPLICIT  REAL*8  (A-H.O-Z) 

DIMEISIOI  A(4) ,B(4) ,C(4) ,XR(3) ,XI(3) 

Ra-10D0 

S=-10D0 


A-4 


MAXIT=100 

S  IF(I.GE.3. AID.ITER.LT.MAXIT)THEI 
ITER=0 

DO  IS  J=1,HAXIT 
ITER=ITER+1 
B(I)=A(I) 

B(I-1)=A(I-1)+R*B(I) 

C(I)=B(I) 

C(I-1)=B(I-1)+R*C(I) 

DO  10  1=1-2, 0,-1 

B(I)=A(I)+R*B(I+l)+S*B(I+2) 

10  C(I)=B(I)+R*C(I+1)+S*C(I*2) 

DET=(DABS(C(2)))**2D0-C(3)*C(1) 

IF(DET. IE. ODO)THEI 
DR=(-n(l)*C(2)+B(0)*C(3))/DET 
DS=(-B(0)*C(2)+B(l)*C(l))/DET 
EA1=DABS(DR/R)*100D0 
EA2=DABS(DS/S) *10000 
R=R+DR 
S=S+DS 
ELSE 
R=R+1D0 
S=S+1D0 
ITER=0 
ERDIF 

IF(EA1.LT. Id-02 . AID . EA2 . LT . Id-02) GOTO  17 
15  COITIIUE 

17  CALL  ROOTS(I,R,S,XR,XI,IEG) 

I=*-2 

DO  20  11=0 ,1 
20  A(II)=B(II+2) 

EIDIF 

IF(I.Gb.3)G0T0  S 
IF ( ITER . LT . MAXIT) THE! 

IF(I.EQ.2)THEI 

R=A(1)/A(2) 

S=A(0)/A(2) 

CALL  ROOTS (I , R , S , XR , XI , IEG) 

ELSE 

XR(1)=-A(0)/A( 1) 

EIDIF 

ELSE 

VRITB(*,30) 

30  FORNAT(SX, 'MAXIMUM  ITERATIOIS  EXCEEDED ! • ,/) 

EIDIF 

RETURI 

EID 

C 

C 

c  M*MM**MM***»*MM»*»*MMMMM*MMM**MMMM* 

C  *  SUBROUTIIE  ROOTS  • 
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c  • 

C  *  THIS  SUBROUTINE  CALCULATES  THE  ROOTS  OF  THE 

C  *  QUADRATIC  USING  R  AND  S  VALUES  CALCULATED  BY 

C  *  SUBROUTINE  POLY. 

c 

SUBROUTINE  ROOTS(N,R,S,XR,XI,NEG) 

IMPLICIT  REAL*8  (A-H.O-Z) 

DIMENSION  XR(3),XI(3) 

F=(DABS(R))**2DO+4*S 

NEG=0 

IF(F.LT.ODO)THEN 

F=-F 

NEG-1 

ENDIF 

XR(H)s(R+DSQRT(F))/2D0 
XR(N-l) * (R-DSQRT(F) ) /2D0 
IF(NEG.EQ. 1)THEN 
XR(N)-R/2D0 
XR(I-1)=R/2D0 
XI(N)=DSQRT(F)/2D0 
XI ( N- 1 ) =-DSQRT(F)/2DO 
ENDIF 
RETURN 
END 
C 
C 

C  ***************************************************** 

c  *  SUBROUTINE  OBSHOC 

C  * 

C  *  THIS  SUBROUTINE  CALCULATES  THE  MACH  NUMBER  ON 

C  *  THE  OTHER  SIDE  OF  A  OBLIQUE  SHQCX  WAVE  AND  THE 

C  *  STATIC  PRESSURE  RATIO. 

C  ***************************************************** 

c 

SUBROUTINE  OBSHOC  (S.MX.MY.PSYX.PTYX) 

IMPLICIT  REAL*8  (A-Z) 

COMMON  K, RAD, DEG 

MY1»(I+1D0)**2D0 

MT2«MX**4D0 

MY3»MX**2DO 

MY4«DSIN(S)**2D0 

MY5-K-1DO 

MY8«(MT1*ICT2*MY4)-4D0*(MY3*MY4-1D0)*(K*MY3*MY4+1D0) 

NY7«(2D0*K*HY3*HY4-NY5)*(MYS*NY3*MY4+2D0) 

MY»DSQRT(MY6/MT7) 

PSYX» (2D0*K*MY3*HY4-MYS ) /(K+ IDO ) 

M1X*HX*DSIN(S) 

PZ1«( (K+ID0)/2D0)*M1X**2D0 
PZ2«l+( (K-1D0)/2D0)*M1X**2D0 
PZ3«IC/(JC-1D0) 
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PZ4=1D0/(1D0-K) 

PZ6s(K-lD0)/(K+lD0) 

PZ6»2D0*K*M1X**2D0/(K+1D0) 

PTYX=((PZ1/PZ2)**PZ3)*(PZ6-PZ5)**PZ4 

RETURI 

EVD 

C 

C 

SUBROUTIIE  SAMAX(M.SAM) 

IMPLICIT  REAL*8  (A-Z) 

COMMOI  K, RAD, DEG 

SAM1=3D0*(3D0*M**4D0+4D0*M**2D0+20D0) 

SAM= ( IDO/ ( 7D0*M**2D0) ) * ( 3D0*M**2D0-SD0+DSQRT (SAMI ) ) 
SAM=DSQRT(SAM) 

SAM=D ASII(SAM) 

RETURI 

EID 

C 

C 

SUBROUTIIE  UAMAX(SA,M,WAM) 

IMPLICIT  REAL*8  (A-Z) 

COMMOI  K.RA.DEG 

WAM1=M**2D0*( (DSII(SA) ) **2D0)-iD0 
WAM=DTAI(SA)*(((6DO*M**2DO)/(SDO*(WAM1)))-1DO) 

UAM=D ATAI ( 1DO/WAM) 

RETURI 

BID 
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A. 2  Data  Aquisition  Program 

The  following  GW  Basic  program  was  written  to  control  the  aquisition  of  the 
pressure  data  through  the  Qua  Tech  Simultaneous  Sampling  System.  The 
Qua  Tech  “LABSIM”  (11)  program  is  listed  in  the  begining  of  the  program  to 
establish  continuity  between  the  Qua  Tech  system  and  the  CPU.  The  remainder 
of  the  program  specially  developed  to  handle  the  unique  data  aquisition  system 
configuration. 


>  * 
>  * 


*  * 
>  * 


LASER  CAVITY  FLOW  DATA  AQUISTIO*  PROGRAM  (LCFDAP) 
WRITTE*  BY:  CAPT  STEPHEI  W.  STIGLICH,  JR. 


WRITTEI  FOR:  THESIS  REQUIREMEITS 
THESIS  ADVISOR:  DR.  V.  C.  ELROD 

This  program  gather*  data  using  a  Qua  Tech  board  thru  pressure 
transducers  located  in  both  stilling  chambers  and  Test  Section. 


**********•*******«*********«**********************»**+•***+********»* 


********* ******************************* 


5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101  ’  ********************************* 

102  » 

220  ADC . SETUP*AH3 • SETCTM-RH6 
29'  SBTC0**H9 

IIADC12.B*AHC:SEGADDR=AHF :  II ADC  12.  BB=AH12 
250  SETMD1**H15 : SETMD2=AH18: SETC1=AH1B : SETC2=AH1E 
260  RSETC1=*H21 :  RSETC2=AH24 : RE ADC 1 =AH27 : READC2=AH2A 
270  DEF  SEG=0 

280  CSEGl*PEEX(RH3CA)*2S6*PEEX(ftH3CB) 

380 
390 
400 

410 

411 

412 


qUA  TECH  ASSEMBLY  LAIGUAGE  R0UTIIES  • 

* 

THE  FIRST  400  LUES  ARE  USED  TO  * 
DEFIVE  AID  LOCATE  SPECIAL  ASSEMBLY  * 
LAIGUAGE  ROUTIIES  REQUIRED  TO  RUI  * 
THE  QUA  TECH  BOARDS.  SEE  THE  QUA  * 
MAIUALS  FOR  A  DESCRIPTIOI  OF  THE  * 
THE  VARIABLES.  * 


* 


APPLICATION  PROGRAM 


420  * 

430  ’  Before  the  CALL  Statements  from  each  program  you  must 
440  *  redefine  the  "SEG"  statements. 

460  ’  For  Labsim  statements  use  DEF  SEG=CSEG1. 

480  * 

490  '************************************************************* 

500  ’  LASER  CAVITY  SET-UP  Six  channel  simultaneous  sampling. 

510  '  The  result  of  the  data  sure  stored  in  the  array  sequentially 

520  ’  array‘/,(l)  —  channel  0 

530  ’  array*/,(2)  —  channel  1 

540  '  arrayX(3)  —  channel  2 

550  '  array*/, (4)  —  channel  3 

560  '  arrayy,(5)  —  channel  4 

570  ’  array!4(6)  —  channel  5 

580  ' 

640  ’ ************************************************************* 

1000  ' 

1005  DIM  CHZV#(6),CHV#(6),P0SITI0IX(6),SL0PEPT#(11) 

1010  COLOR  7,1: CLS  0 

1020  PRIIT  "  DATA  RUI  COIFIGURATIOI" 

1025  PRIIT  "  MEIU" 

1030  PRIIT  "" 

1035  PRIIT  "COIFIG  #  TRAISDUCER  POS.  #  ASSIGIED  TO  QUA  TECH  CHAIIEL 

#" 


1040 

1045 

PRIIT 

PRIIT 

M 

MH 

0 

1 

2 

3 

4 

5" 

1050 

PRIIT 

II 

1 

01 

02 

03 

04 

05 

06" 

1055 

PRIIT 

II 

2 

01 

02 

15 

12 

05 

06" 

1060 

PRIIT 

M 

3 

01 

17 

14 

11 

09 

06" 

1065 

PRIIT 

II 

4 

18 

16 

03 

13 

10 

06" 

1070 

PRIIT 

II 

5 

07 

14 

15 

16 

08 

06" 

1075 

PRIIT 

II 

6 

17 

02 

18 

13 

11 

06" 

1080  PRIIT  "  7  OPERATOR  ASSIGIED." 

1086  PRIIT  "" 

1090  IIPUT  "SELECT  THE  COIFIGURATIOI  FOR  THIS  RUI.  ".C0IF7. 

1100  IF  C0IF%*1  THEI  GOSUB  10000 
1110  IF  C0IF%=2  THEI  GOSUB  20000 
1120  IF  C0IFX*3  THEI  GOSUB  30000 
1130  IF  COIFX-4  THEI  GOSUB  40000 
1140  IF  C0IF%*6  THEI  GOSUB  50000 
1150  IF  C0IF%*6  THU  GOSUB  60000 
1160  IF  COIFX-7  THEI  GOSUB  65000 
2000  * 

2010  CLS  0 

2015  SL0PEPTi(2)*2. 39741157912 

2016  SL0PEPT#(8)«3. 11210203183 

2020  SAMPLERATE%*200 

2021  SEC* IS 

2030  IIPUT  "DO  YOU  WAIT  TO  TAKE  A  IEW  ZERO  RUI  (YES*Y,  10=1) ";Z$: PRIIT  "" 

2031  IF  Z8«"I"  OR  Zl*"n"  THEI  GOTO  2100 

2032  IIPUT  "EITER  FILEIAME  FOR  ZERO  RUI  DATA:  ",  ZDlt: PRIIT 


II II 


2033  ZDl!="C : \SWS\DATA1\"+ZD l!+" . DAT" 

2080  COLOR  19, 1  :CLS  0:PRIIT  "":PRIIT  ,,M:PRIIT  ""-.PRIIT  •••• 


2070  PRIIT  M  ****•*************•**********♦♦♦*,*» 

2071  PRIIT  "  *  ZERO  VOLTAGE  RUI  II  PROGRESS!  *" 

2072  PRIIT  "  ****************************»******m 


2073  ZDATEt=DATE! 

2074  ZTIME!=TIME! 

2075  GOSUB  5000 
2080  ’ 

2081  OPEI  ZDl!  FOR  OUTPUT  AS  #1 

2082  PRIIT  #1."  ZERO  VOLTAGE  RUI*' 

2083  PRIIT  #1,M" 

2084  PRIIT  #1,'‘DATE  OF  RUI:  “;ZDATE! 

2085  PRIIT  #1 ."TIME  OF  RUI:  ";ZTIME! 

2086  PRIIT  #1,"" 

2087  FOR  1=1  TO  6 

2088  PRIIT  #1 ,CHZV#(I) 

2089  I EXT  I: CLOSE  *1 

2090  COLOR  31 , 1  :CLS  0: PRIIT  ""-.PRIIT  "": PRIIT  "": PRIIT 


2091  PRIIT  "  ***********************************>• 

2092  PRIIT  ”  *  ZERO  VOLTAGE  RUI  COMPLETED!  *" 

2093  PRIIT  "  ************************************* 

2094  PRIIT 


2095  PRIIT: COLOR  7.1 

2096  GOTO  2110 

2100  IIPUT  "EITER  FILEIAME  OF  ZERO  DATA  FILE  TO  USE:  ",ZD2!: PRIIT  — 

2101  ZD2l="C:\S¥S\DATAl\»+ZD2t*". DAT- 

2102  OPEI  ZD2$  FOR  IIPUT  AS  «1 

2103  FOR  1=1  TO  5 

2104  LIIE  IIPUT  #1,  JUIKI 

2105  I EXT  I 

2106  FOR  1=1  TO  6 

2107  IIPUT  #1.  CHZVi(I) 

2108  IEXT  I 

2109  CLOSE  *l:ZDlt=ZD2I 

2110  IIPUT  -FILEIAME  FOR  DATA  (EXTEISIOI  HILL  BE  .DAT)"; It: PRIIT  — 

2113  It«"C : \SWS\D ATA 1 \" ♦!!♦" .DAT- 

2114  IIPUT  "EITER  THE  DIFFUSER  THROAT  WIDTH  (IICHES)" ;TW#:PRIIT  "" 

2115  PRIIT  "TEST  LEI GTE  ■  ";SEC;"  SECOIDS" : PRIIT  — 

2116  IIPUT  "DO  YOU  WAIT  TO  CHAIGE  THE  TEST  LEIGTH  (YES»Y,  10=1) -;TLt: PRIIT 

M  II 

2117  IP  TU--I"  OR  TLI»"n"  THE!  2119 

2111  IIPUT  "BITER  KW  TEST  LEIGTI.  ",  SEC: PRIIT  — 

21../  PRIIT  "SAMPLE  UTE  ■  " ;SAMPL£RATE%: PRIIT  — 

2120  IIPUT  "DO  YOU  WAIT  TO  CHAIGE  THE  SAMPLB  RATE  (YES«Y,  10=1) ";SRt: PRIIT 

M  II 

2121  IF  SRt»"I"  OR  SRt»"n"  THE!  2125 

2122  IIPUT  "EITER  IEW  SAMPLE  RATE.  ", SAMPLERATEX : PRIIT  — 

2126  IIPUT  "TURI  01  LAB  VACUUM  PUMPS,  THE!  PRESS  RETURI. »; DUMMY!: PRIIT  — 

2126  IIPUT  "TURI  01  REP  VACUUM  PUMP,  THEI  PRESS  RETURI DUMMY! : PRIIT  — 

2127  IIPUT  -EITER  REPEREICE  PRESSURE  (MM  OP  HG):  PREF# : PRIIT  — 
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2128  PREF#=PREF#*. 019264 
2130  GOSUB  8000 

2150  IIPUT  "EITER  ATMOSPHERIC  PRESSURE  (II  OF  HG) :  H;PATM#:PRI*T 
2180  IIPUT  "EITER  REFEREICE  PRESSURE  (MM  OF  HG) :  " ; PREF# : PRUT  "" 

2170  PATM#=PATM#* . 4893 
2180  PREF#=PREF#*. 019264 

2190  POIITSX=SEC*SAMPLERATEX 

2191  DIMAX=P0IITSX*6 

2192  TESTDATE$=DATE$ 

2193  IIPUT  "PRESS  RETURI  TO 

2195  TESTTIME$=TIME| 

2196  COLOR  19,1:CLS  0:PRIIT 

2197  PRUT  " 

2198  PRUT  " 

2199  PRUT  " 

2200  GOSUB  6000 

2201  COLOR  31 , 1 :CLS  0:PRIIT 

2202  PRUT  " 

2203  PRUT  " 

2204  PRUT  " 

2206  PRUT  " 

2207  PRUT  " 

2208  PRUT 

2209  PRUT: COLOR  7,1 

2210  OPEI  1$  FOR  OUTPUT  AS  #1 

2211  PRUT  #1."  TEST  DATA  RUI" 

2212  PRUT 

#1 ,P0SITI0IX(1) ,P0SITI0iy.(2) ,P0SITI0I%(3) ,P0SITI0IX(4) ,P0SITI0I%(5) , POSITIO 
IX(6) 

2213  PRUT  #1, SEC, TV# 

2216  PRIIT  #1 , SAHPLERATEX 

2217  PRIIT  *1 .TESTDATEl .TESTTIMEl 

2218  PRIIT  #1 ,CHZV#(1)  CHZV#(2)  CHZV#(3)  CHZV#(4)  CHZV#(5)  CHZV#(6) 

2219  PRIIT  # 1 , P ATM# , PREF# 

2220  IIDEXX*1 

2221  FOR  J=1  TO  DIMAX 

2222  PRIIT  #1.ARRAYX(J);"  "; 

2223  IIDEX%*IIDEXX+1 

2224  IF  IIDIXX«7  THEI  PRIIT  #1,"" 

2226  IF  IIDEXX-7  THEI  IIDEXX=1 

2228  IBXT  J 

2229  CLOSB  #1 

2300  HISZ$*"C: \SVS\HIS\HISZ.DAT" 

2310  0PS1  HISZ8  FOR  APPEID  AS  #1 
2320  PRIIT  #1,ZD19 
2330  CLOSE  #1 

2400  HISTt*"C : \SVS\HIS\HIST.DAT" 

2410  OPEI  HIST8  FOR  APPEID  AS  #1 
2420  PRIIT  #1,1$ 

2430  CLOSE  #1 
2900  CLS  0 


START  DATA  RUI.";DUMMY8 

PRIIT  "": PRIIT  "": PRIIT  "" 
***********************************" 

*  DATA  RUI  II  PROGRESS!  *" 
***********************************" 

"": PRIIT  ""-.PRIIT  "": PRIIT  "" 

***********************************" 

*  DATA  RUI  COMPLETE!  *" 

*  »*• 

*  PRIITIIG  DATA  TO  FILE  * 

***********♦***********************" 
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2920  COLOR  15,1:  PRUT  m,:PRIHT  MM 
2930  SOUID  2600,10 

2950  IIPUT  "DO  YOU  WAIT  TO  TAKE  AIOTHER  DATA  RUI?  " ; AIS* : CLS  0 

2960  IF  AIS$-"IH  OR  AIS$="n"  THE!  GOTO  3000 

2970  ERASE  ARRAYS: GOTO  1010 

3000  CLS  0: CLEAR: SYSTEM 

5000  ' 

5010  '«*************#************#****************•***•••*•*•*****« 
5020  >*  SUBROUTIIE  ZERO  * 

5030  '*  * 

5040  '*  This  subroutine  samples  the  pressure  transducers  and  * 

5050  ’*  calculates  the  zero  voltages  prior  to  taking  a  run.  * 

5060  > tee********************************************************** 
5070  ' 

5100  DEF  SEG=CSEG1  ’DEFIIE  THE  A  TO  D  CARD 

5110  DIM  ARRAYX(20000) 

5120  ADDRESSX=*H300 

5130  I7,=6  '#  OF  CHAIIELS 

5140  CALL  ADC .  SETUP  (  ADDRESS'/, ,  1%) 

5150  C0%=5000000 ! /20000  ' 20KHZ  SAMPLIIG  RATE 

5160  CALL  SETCO(COX) 

5170  LEIGTHX^lOO  'TAKE  100  POUTS 

5180  LEIGTH%=LEIGTHX*6 
5190  S#=LEIGTH%/6 

5200  CALL  IIADC12 . B (LEIGTH7, , ARRAYXC 1 ) )  'TAKE  ZERO  VOLTAGES  DATA 

5210  FOR  1=1  TO  6 
5220  CHZV#(I)=0! 

5230  IEXT  I 
5270  FOR  J=1  TO  6 

5276  FOR  I=1+(J-1)  TO  S9S+(J-1)  STEP  6 

5280  CHZV#(J)=(2047-ARRAYX(I))/819.2+CHZV#(J)  'COIVERT  TO 

VOLTAGES 

5340  IEXT  I 

5346  IEXT  J 

5350  ' 

5360  ’ 

5370  ' 

5380  FOR  1=1  TO  6 

6390  CHZYt(I)  •  CHZV#(I)/S#  'PSIG  TRAISDUCER 

5400  IEXT  I 
5436  * 

5460  ERASE  ARRAYX 
5466  RSTURI 
6000 
8010 
6020 
6030 
6040 
6050 
6060 
6070 


e  SUBROUTIIE  PRES 

• 

•  This  subroutine  saaples  the  pressure  transducers  and 

*  calculates  the  pressure  voltages. 
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'DEFIIE  THE  A  TO  D  CARD 


6100  DEF  SEG-CSEG1 
6110  DIN  ARRAYy.(DINAX) 

6120  ADDRESSX=AH300 
6130  IX= 6  ’#  OF  CHA1IELS 

6140  CALL  ADC. SETUP (ADDRESSy.,My.) 

6150  C0y.=5000000!/SANPLERATEy.  '  5MHZ/SAMPLIIG  RATE 

6160  CALL  SETC0(C0y.) 

6170  LEIGTHy.=POI*TSy.  'TAKE  SO  MAIY  POUTS  PER  CHHL 

6200  CALL  I*ADC12.B(LE*GTH*/.,ARRAYX(1))  'TAKE  PRES  VOLTAGES  DATA 

6470  RETURI 
8000  ' 

8001  '****************  ******************************************** 

8002  '♦  SUBROUTIIE  TEST  SECTIOI  PRESSURE  * 

8003  '*  * 

8004  ' *  This  subroutine  samples  the  pressure  transducer  in  the  * 

8005  **  test  section  and  gives  the  operator  a  pressure  value.  * 

8006  '*  This  allows  the  operator  to  monitor  the  test  section  * 

8007  **  pressure  and  determine  when  he  is  ready  to  make  a  run.  * 

8008  ************************************************************** 

8009  ' 

8025  CLS  0 

8060  VIEW  PRIIT  1  TO  IS 
8055  SL0PE#=SL0PEPT#(2) 

8060  IF  P0SITI0*y.(2)  >  2  THE*  SL0PE#=SL0PEPT#(8) 

8080  LOCATE  6,20 
8090  COLOR  16,1 

8100  PRIIT  "*************************************" 

8105  LOCATE  6,20 

8110  PRIIT  »*  PRESS  RETURI  OICE  VACUUM  I*  TEST  *M 
8115  LOCATE  7,20 

8120  PRIIT  "*  SECTIOI  HAS  STABILIZED.  *" 

8125  LOCATE  8,20 

8130  PRIIT  ”*************************************" 

8200  VIEW  PRIIT  11  TO  12 
8210  COLOR  15,1 

8300  P$=IIKEYt:IF  LE*(P$)<>0  THE1  8900 

8310  DEF  SEG=CSEG1  'DEFIIE  THE  A  TO  D  CARD 

8320  DIN  ARRAYXC 20000) 

8330  ADDRESS%**H300 

8340  IX«6  '*  OF  CHAIIELS 

8360  CALL  ADC. SETUP (ADDRESSX, *%) 

8360  COX«6000000!/20000  ' 20KHZ  SAMPLING  RATE 

8370  CALL  SETCO(COX) 

8380  LEIGTHX>100  ‘TAKE  100  POIITS 

8390  LEIGTHX-LEIGTHX*IX 
8400  S*>LEIGTHX/IX 

8410  CALL  I*ADC12.B(LE*GTHX,ARRAYX(1))  'TAKE  ZERO  VOLTAGES  DATA 

8420  FOR  1-1  TO  IX 
8430  CHV#(I)*0! 

8440  I EXT  I 

8450  FOR  J*1  TO  IX 
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8480  FOR  I=1+(J-1)  TO  595+(J-i)  STEP  6 

8470  CHV#(J)=(2047-ARRAYy.(I))/819.2!+CHV#(J)  'COIVERT  TO 

VOLTAGES 

8480  I EXT  I 

8490  I EXT  J 

8500  ' 

8510  ' 

8520  » 

8530  FOR  1=1  TO  6 

8640  CHV#(I)  =  CHV#(I)/S#  >PSIG  TRANSDUCER 

8550  I EXT  I 
8560  ' 

8570  ERASE  ARRAYS 

8600  CHV#(2)=CHV#(2)-CHZV#(2) 

8610  PCH2G#=CHV# ( 2 ) *SLOPE# 

8620  PCH2A#=PREF#+PCH2G# 

8710  PRUT  "TEST  CHAMBER  STAGIATIOV  PRESSURE  IS: 

8715  PRIIT  USIIG  "#***.***#«" ;PCH2A«; 

8716  PRIIT  "  PSIA" 

8718  IF  PCH2A#  <  .1350  THEI  SOUVD  2600,10 
8720  GOTO  8300 

8900  VIEW  PRIIT  1  TO  26:CLS  0 
8950  RETURI 
10000  * 

10010  P0SITI0lX(l)=01 
10020  P0SITI0I%(2)-02 
10030  P0SITI0I% ( 3 ) =03 
10040  P0SITI0IX(4)=04 
10050  P0SITI0I%(5) =05 
10060  P0SITI0I%(6)*06 
10100  RETURI 
20000  > 

20010  P0SITI0IX(1)=01 
20020  P0SITI0I%(2)*02 
20030  P0SITI0IX(3)>16 
20040  P0SITI0I%(4)«12 
20050  P0SITI0I%(5)*06 
20060  P0SITI0IX(6)«06 
20100  RETURI 
30000  • 

30010  P0SITI0IX(1)>01 
30020  P0SITI0IX(2)«17 
30030  P0SITI0I%(3)«14 
30040  P0SITX0IX(4)«11 
30060  P0SITI0VX(6)«09 
30060  P0SITI0I%(6)*06 
30100  RETURI 
40000  * 

40010  P0SITI0IX(1)>18 
40020  P0SITX0K(2)*16 
40030  P0SITI0I%(3)*03 
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40040  P0SITI0IX(4)=13 
40050  P0SITI0I%(5)=10 
40060  POSITIOI* (6) =06 
40100  RETURI 
50000  » 

50010  P0SITI0I%(1)=07 
50020  P0SITI0I*(2)=14 
50030  P0SITI0I%(3)=15 
50040  P0SITI0IX(4)=16 
50050  P0SITI0I%(5)=08 
50060  P0SITI0I%(6)=06 
50100  RETURI 
60000  » 

60010  P0SITI0I%(1)=17 
60020  P0SITI0I*(2)=02 
60030  P0SITI0iy.(3)=18 
60040  P0SITI0«y.(4)=13 
60050  POSITIO*y.(5)  =  ll 
60060  P0SITI0I%(6)=06 
60100  RETURI 
65000  ' 

65002  CLS  0 


65010  PRUT 
65020  PRUT 
66030  PRUT 
65040  PRUT 
85050  PRUT 
65060  PRUT 
65070  PRUT 
65075  PRIIT 
66076  PRIIT 

65077  PRIIT 

65078  PRIIT 
66079  PRIIT 
65080  PRIIT 
65085  PRIIT 
65090  PRIIT 
65096  IIPUT 


OPERATOR  ASSIGIED  COIFIGURATIOI" 
MEIU" 

CURREIT  QUA  TECH  CHIL  ASSIGIMEITS: 


CHAIIEL  # 


POSITIOI  #" 

POSITIOIX(I) 

P0SITI0iy.(2) 

P0SITI0I%(3) 

P0SITI0iy.(4) 

POSITIOIX(S) 

P0SITI0I%(6) 


65090  PRIIT  HEITER  CHAIIEL  *  TO  CHAIGE  OR  EITER  9  VHEI  COMPLETE." 

65096  IIPUT  "REMEMBER  TO  IIPUT  TWO  DIGITS  (EXAMPLE:  01  OR  18).  ",OPC% 
66100  IF  0PC%»0  TIBI  IIPUT  "EITER  POSITIOI  ASSIGIMEIT  FOR  CHAIIEL  0. 
".POSITIQtt(l) 

66110  IF  OPCX-1  THBI  IIPUT  "EITER  POSITIOI  ASSIGIMEIT  FOR  CHAIIEL  1. 

" ,P0SITI0I%(2) 

66120  IF  0PCX>2  THU  IIPUT  "EITER  POSITIOI  ASSIGIMEIT  FOR  CHAIIEL  2. 
",P0SITI0rC(3) 

65130  IF  0PCX>3  THEI  IIPUT  "EITER  POSITIOI  ASSIGIMEIT  FOR  CHAIIEL  3. 

" ,P0SITI0I%(4) 

65140  IF  0PC%*4  THEI  IIPUT  "EITER  POSITIOI  ASSIGIMEIT  FOR  CHAIIEL  4. 

" ,P0SITI0lX(6) 

65160  IF  0PC%>5  THEI  IIPUT  "EITER  POSITIOI  ASSIGIMEIT  FOR  CHAIIEL  5. 

" ,P0SITI0I%(6) 

65160  IF  0PC%>9  THEI  GOTO  65180 
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65170  GOTO  65000 
65180  RETURI 
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A. 3  Data  Reduction  and  System  Control  Program 

This  Quick  Basic  program  controls  of  both  the  data  aquisition  and  the  data 
reduction.  It  also  provides  the  capability  to  plot  the  reduced  data  files.  The 
operator  needs  only  enter  Laser  and  follow  the  procedures  outline  on  the  screen. 
Subroutine  Axis  was  provided  by  Tanis  (15). 


DECLARE  SUB  VOLTAGE  (CH!(),  ARRAY*/, () ,  DIMAX,  IX.  FLAGX,  TX.  TSX,  IlX,  I2X) 
DECLARE  SUB  PLOTMEVU  (TDC$,  CHILI) 

DECLARE  SUB  CHA1IEL  (TDCl,  CHILI.  TIME!(),  PRES!()) 

DECLARE  SUB  MAIIMEIU  () 

DECLARE  SUB  DATARUI  () 

DECLARE  SUB  DATACOIV  () 

DECLARE  SUB  PLOTDATA  () 

DECLARE  SUB  MAXIMUM  (max,  DMAX!(),  TSTEPX) 

DECLARE  SUB  AXIS  (XMAX!,  XMIV!,  YMAX!  ,  YMII!) 

DECLARE  SUB  MIIIMUM  (■in,  DMII!(),  TSTEPX) 

COMMOI  SHARED  SAMPLERATEX,  IX,  SEC,  TSTEPX 

SCREE!  9 

COLOR  7,  1 

SEC  =  15 

CALL  MAIIMEIU 

SCREE!  0 

CLS  0 

COLOR  1,  1 
EID 

SUB  AXIS  (XMAX,  XMII,  YMAX,  YMII)  STATIC 
DIM  X(0  TO  8),  y(0  TO  8) 

EGAl  *  "Y" 

VIEW  PRIIT  1  TO  24 
DX  =  (XMAX  -  XMII)  /  5 
DY  »  (YMAX  -  YMII)  /  5 
FOR  I  a  0  TO  5 

X(I)  ■  I  *  DX  +  XMII 
y(I)  »  I  *  DY  +  YMII 

■EXT  I 

FOR  I  »  5  TO  0  STEP  -1 

LOCATE  18  -  I  *  3,  1 
PRIIT  USIIG  . . ;  y(I) 

■EXT  I 
XI  «  7 

LOCATE  20,  1 
PRIIT  TAB(XI) ; 

FOR  I  a  0  TO  5 

PRIIT  TAB (XI  ♦  13  •  I); 

PRIIT  USIIG  "#.## -  X(I) ; 
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I  EXT  I 
PRUT 

'SETS  VIEW  PORT  FOR  CGI 

IF  EGA!  =  "Y"  THEM  VIEW  (70,  3S)-(S90,  254)  ELSE  VIEW  (70, 
20)-(S90,  145) 

'DRAWS  TIC  MARKS  FOR  THE  AXIS 
IF  EGAS  *  "YM  THE! 

DRAW  "BH77,35;IL7" 

FOR  I  =  0  TO  4 

DRAW  "BM+0,11;IL3" 

DRAW  "BM+O, 10;IL3" 

DRAW  "BM+O, 11;IL3" 

DRAW  "BM+O, 10;IL7" 

I  EXT  I 

ELSE  DRAW  "BM77,20;IL7" 

FOR  I  =  0  TO  4 

DRAW  "BM+0,6;IL3" 

DRAW  "BM+0,6;*L3" 

DRAW  "BM+0,6;IL3" 

DRAW  "8M+0,fl;*L7" 

I  EXT  I 
EVD  IF 

IF  EGAS  a  "Y"  THE*  DRAW  »BM80,247;ID6"  ELSE  DRAW  »BM80, 141 ; *D8" 
FOR  I  =  0  TO  5 

DRAW  "BM+25,0;ID2" 

DRAW  "BM+26 , 0 ; *D2" 

DRAW  "BM+25,0;*D2" 

DRAW  "BM+26, 0;*D7" 

IEXT  I 

'SETS  VIEW  PORT  FOR  PLOTIIG  THE  DATA 

IF  EGAS  »  "YM  THE* 

VIEW  (80,  35)-(590,  245),  9,  11 
ELSE  VIEW  (80,  20)-(590,  140),  ,  2 
EID  IF 

'SETS  SCALE  OF  THE  VIEW  PORI 

'WIIDOW  (-2,0)-(6.5,4.7S) 

WIIDOW  (XMX1,  YMX1)-(XMAX ,  YMAX) 

'  LIIE  (XMII ,0)-(XMAX ,0) 

VIEW  PRIIT  21  TO  24 

EID  SUB 

SUB  CHAIIEL  (TDCS,  CHILS,  TIME'O,  PRES!()) 

plilS  ■  "C:\SWS\DATA1\"  +  TDCS  +  CHILS  +  ".DAT" 

OPEI  plilS  FOR  IIPUT  AS  *20 
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IIPUT  #20,  TSTEP‘/. 

FOR  I  =  1  TO  TSTEPX 

IIPUT  #20,  TIME! (I) 

IIPUT  #20,  PRES! (I) 

■EXT  I 
CLOSE  #20 
EID  SUB 

SUB  DATACOIV 

DIM  CH ! (8) ,  CHZV! (6) ,  PCHG!(6),  PCHA!(6),  ARRAY* (1 8000 ) ,  SL0PEPT!(11), 
P0SITI0iy.(8) 

CLS  0 

IIPUT  "Ester  the  number  o 2  samples  per  point  on  the  graph.  ",  TS* 
HISS  =  "C:\SVS\HIS\HIST.DAT" 

iy.  =  e 

SCREEI  0 
COLOR  7,  1 
CLS  0 

OPEI  HISS  FOR  IIPUT  AS  #1 
DO  UITIL  EOF(l) 

LI IE  IIPUT  #1,  DATAFILIAMS 
FILIAMES  «  RIGHTS (DATAFILIAMS,  10) 

CLS  0 
COLOR  23 
LOCATE  8.  25 

PRUT  “READIIG  DATA  FILE  "  +  FILIAMES  + 

TESTDS  *  "C:\SVS\DATA1\"  ♦  FILIAMES  +  "" 

OPEI  TESTDS  FOR  IIPUT  AS  #2 
FOR  I  *  1  TO  7 
IF  I  a  1  THE! 

LIRE  IIPUT  #2,  JUIKS 
EID  IF 

IF  I  =  2  THE! 

IIPUT  #2,  POSITIOIX(l),  P0SITI0lX(2) ,  P0SITI0*X(3) , 
P0SITI0I*(4) ,  POSITIOIX(S) ,  P0SITI0I%(8) 

EID  IF 

IF  I  a  3  THE! 

IIPUT  #2,  SEC,  TV! 

EID  IF 

IF  I  «  4  THE! 

IIPUT  #2,  SAMPLERATEX 
EID  IF 

IF  I  a  5  THEI 

IIPUT  #2,  DATETIMES 
DATES  «  LEFTS (DATETIMES,  10) 

TIMES  «  RIGHTS (DATETIMES.  8) 

EID  IF 

IF  I  a  6  THEI 

IIPUT  #2,  CHZV ! ( 1 ) ,  CHZV! (2) ,  CHZV!(3),  CHZV!(4),  CHZV!(5), 

CHZV! (8) 

EID  IF 
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IF  I  =  7  THE* 

IIPUT  #2,  PATH! ,  PREF! 

EID  IF 
IEZT  I 

POIITSX  =  SEC  *  SAMPLERATE* 

DIMA*  =  POUTS*  *  I* 

FOR  I  =  1  TO  DIMA* 

IIPUT  #2,  ARRAY*(I) 

IEXT  I 
CLOSE  *2 
CLS  0 
COLOR  23 
LOCATE  8,  25 

PRIIT  "COIVERTIIG/SAVIIG  DATA  FILE  "  +  FILIAMES  + 

SLOPEPT! (1)  =  3.532384 

SLOPEPT! (2)  =  2.397412 

SLOPEPT! (3)  =  3.333312 

SLOPEPT! (4)  =  2.804098 

SLOPEPT! (5)  =  12.81 

SLOPEPT! (6)  -  63.58949 

SLOPEPT! (7)  =  4.205204 

SLOPEPT! (8)  *  3.113102 

SLOPEPT! (9)  *  3.069795 

SLOPEPT! (10)  *  2.804098 

SLOPEPT! (11)  *  3.872598 

TSTEPX  =  (SAMPLERATEX  *  SEC)  /  TS* 

TD$  -  LEFTS (FILIAMES ,  8) 

P0ST1S  =  LTRIMS ( STRS ( POSITIOI* ( 1 ) ) ) 

P0ST2S  *  LTRIMS (STR$(P0SITI0I*(2))) 

P0ST3S  *  LTRIMS ( STRS ( POSITIOI* ( 3 ) ) ) 

P0ST4S  =  LTRIMS (STRS (P0SITI0**(4))) 

P0ST5S  *  LTRIMS (STRS ( POSITIOI* (5 ) ) ) 

POSTflS  =  LTRIMS ( STRS ( POS ITI 01% ( 6 ) ) ) 

IF  POSITIOIX(I)  <  10  THE* 

P0ST1S  *  "0"  +  P0ST1S  + 

EID  IF 

IF  P0SITI0IX(2)  <  10  THE* 

P0ST2S  •  "0"  ♦  P0ST2S  ♦  t 

EID  IF 

IF  POSITIOI* (3)  <  10  THEI 
P0ST3S  ■  "0"  ♦  P0ST3S  ♦ 

EID  IF 

IF  P0SITI0IX(4)  <  10  THEI 
P0ST4S  ■  ”0"  +  P0ST4S  +  "" 

EID  IF 

IF  POSITIOlX(S)  <  10  THEI 
POSTS*  »  "0"  ♦  P0ST5S  +  "" 

EID  IF 

IF  P0SITI0I*(8)  <  10  THEI 
POSTS*  ■  "0"  +  POSTS*  +  "M 
EID  IF 
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TDCHlI  s  "C:\SWS\DATA1\"  ♦  TDl  ♦  POSTlI  +  ".DAT" 

TDCH2I  =  "C : \SVS\DATA1\"  ♦  TDl  +  P0ST2I  +  ".DAT" 

TDCH3I  =  "C:\SVS\DATA1\"  ♦  TDl  +  P0ST3I  +  ".DAT" 

TDCH4I  *  "C:\SVS\DATA1\"  +  TDl  ♦  P0ST4I  +  ".DAT" 

TDCH5I  »  "C:\SWS\DATA1\"  ♦  TDl  +  P0ST5I  +  ".DAT" 

TDCH6I  =  "C:\SWS\DATA1\"  +  TDl  +  P0ST6I  +  ".DAT" 

OPEI  "C:\SWS\HIS\HISP.DAT"  FOR  APPEID  AS  #3 
PRUT  #3,  ""  +  TDl  +  POST  1 1  +•  »» 

PRIIT  #3,  ""  +  TDl  ♦  P0ST2I  +  "" 

PRUT  #3,  ""  +  TD|  +  P0ST3I  +  "" 

PRIIT  #3,  ""  +  TDl  +  P0ST4I  ♦  "" 

PRIIT  #3,  ""  +  TDl  +  POSTSI  +  "» 

PRIIT  #3,  “"  +  TDl  ♦  P0ST6I  +  »" 

CLOSE  *3 

OPEI  TDCH1I  FOR  OUTPUT  AS  #4 
PRIIT  #4,  TSTEPX 
CLOSE  *4 

OPEI  TDCH2I  FOR  OUTPUT  AS  #4 
PRIIT  #4,  TSTEPX 
CLOSE  #4 

OPEI  TDCH3I  FOR  OUTPUT  AS  #4 
PRIIT  #4,  TSTEPX 
CLOSE  *4 

OPEI  TDCH4I  FOR  OUTPUT  AS  #4 
PRIIT  #4,  TSTEPX 
CLOSE  *4 

OPEI  TDCHSI  FOR  OUTPUT  AS  #4 
PRIIT  #4.  TSTEPX 
CLOSE  *4 

OPEI  TDCHSI  FOR  OUTPUT  AS  94 
PRIIT  #4,  TSTEPX 
CLOSE  *4 

OPEI  "C:\SWS\RPT\"  +  TDl  ♦  ".FIL"  FOR  OUTPUT  AS  #S 

PRIIT  #S,  "  DATA  RUI  »  +  TDl  +  » 

RESULTS" 


PRIIT  #5,  "" 

PRIIT  #6,  "DATE: 

PRIIT  #5,  "TIME: 

PRIIT  95,  "ATMOSPHERIC  PRESSURE: 
P1IIT  95,  " REFERENCE  PRESSURE: 
PRIIT  96,  "DIFFUSER  AREA  RATIO: 
PRIIT  95,  "SAMPLE  RATE: 

PRIIT  95,  "" 

PRIIT  95, 


;  DATE! 

;  TIMEI 
PATM! 

PREF! 

TV! 

SAMPLERATEX 


"TIME  (SEC) 

PRIIT  95,  SPC(IS);  POSTlI;  SPC(7);  P0ST2I;  SPC(7); 
PRIIT  95,  P0ST3I;  SPC(7) ;  P0ST4I;  SPC(7); 

PRIIT  95,  POSTS! ;  SPC(7);  P0ST6I 
PRIIT  96.  "" 

CLOSE  95 

FOR  TX  »  1  TO  TSTEPX 
FLAGX  «  IX  •  TSX 


POSITIOI" 
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CALL  VOLTAGE (CH!  () ,  ARRAYXO,  DIMA*/.,  I*/.,  FLAG’/.,  T X,  TSX,  IlX,  12*/.) 
FOR  I  =  1  TO  *% 

CH!(I)  *  CH !  (I)  -  CHZV!(I) 

I  EXT  I 

FOR  I  =  1  TO  *X 
SLPX  =  I 

IF  POSITIOiy.(I)  >  6  THEM  SLPX  =1  +  6 
PCHG!  (I)  =  CH !  (I)  *  SLOPEPTKSLPy.) 

IEXT  I 

PCHA ! (1)  =  PREF!  +  PCHG!(1) 

PCHA ! (2)  =  PREF!  +  PCHG! (2) 

PCHA! (3)  =  PREF!  +  PCHG! (3) 

PCHA! (4)  =  PATH!  +  PCHG! (4) 

IF  P0SITI0IX(4)  >  4  THE* 

PCHA! (4)  =  PREF!  +  PCHG! (4) 

EID  IF 

PCHA! (6)  =  PATH!  +  PCHG! (6) 

IF  POSITIOIX(S)  >  5  THE* 

PCHA! (6)  =  PREF!  +  PCHG! (5) 

EID  IF 

PCHA! (6)  =  PATH!  +  PCHG! (6) 

TIM!  =  CDBL((1  /  SAMPLERATEX)  *  TX  *  TSX) 

OPE*  TDCHlt  FOR  APPEID  AS  *3 

PRUT  #3,  USIIG  "###.#####";  TIM!; 

PRUT  #3,  USIIG  "###.#####";  PCHA!(1) 

CLOSE  #3 

OPE*  TDCH2I  FOR  APPEKD  AS  *4 

PRUT  #4,  USIIG  "###.#####";  TIM!; 

PRIIT  *4,  USIIG  "###.#####" ;  PCHA! (2) 

CLOSE  *4 

OPEI  TDCH3I  FOR  APPEID  AS  *S 

PRIIT  #5,  USIIG  "###.###••";  TIM!; 

PRIIT  #6,  USIIG  "###.####•";  PCHA! (3) 

CLOSE  *5 

OPEI  TDCH4S  FOR  APPEID  AS  #6 

PRIIT  *6,  USIIG  "**«.*****”;  TIM!; 

PRIIT  *6,  USIIG  "###.#####";  PCHA! (4) 

CLOSE  *6 

OPEI  TDCBSt  FOR  APPEID  AS  97 

PRIIT  f7,  USIIG  "###. #####••;  TIM!; 

PRIIT  91,  USIIG  "###.#####";  PCHA! (5) 

CLOSE  #7 

OPEI  TDCI6$  FOR  APPEID  AS  *8 

PRIIT  #8,  USIIG  "###. #####";  TIM!; 

PRIIT  88,  M  "; 

PRIIT  #8,  USIIG  "#*i. #*##*";  PCHA! (6) 

CLOSE  98 

OPE*  "C:\SWS\RPT\"  +  TDI  ♦  ” . FIL"  FOR  APPEID  AS  #5 
PRIIT  »8.  USIIG  "###. #####";  TIM!; 

PRIIT  #5,  "  "; 

PRIIT  fS,  USIIG  "##i . #####" ;  PCHA ! ( 1) ;  PCHA! (2);  PCHA! (3); 
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PCHA ! (4) ;  PCHA!(5);  PCHJU  (6) 

CLOSE  *5 

PER!  =  CS*G(100  *  tx  /  TSTEPX) 

VIEW  PRIIT  12  TO  13 
CLS 

LOCATE  12,  25 

PRIIT  "File  PER!;  "  X  COMPLETE!" 

VIEW  PRUT  1  TO  25 
■EXT  TX 
LOOP 
CLOSE  *1 

SHELL  "CD\SWS\HIS" 

SHELL  "ERASE  HIST.DAT" 

COLOR  7 
CLS  0 
EID  SUB 

SUB  DATARUI 

SHELL  "CD\SWS" 

SHELL  "LABSIM" 

SHELL  "BASICA  LCF3.BAS" 

CLS  0 
SCREE!  0 
COLOR  7,  1 

EID  SUB 

SUB  MAIIMEIU 
100  CLS  0 

VIEW  PRIIT  1  TO  25 

PRIIT  "  Main  Mann  Selections" 

PRIIT 

PRIIT  "1.  Data  Ron" 

PRIIT  "2.  Data  Conversion" 

PRIIT  "3.  Plot  Data" 

PRIIT  "4.  Exit  Program" 

PRIIT 

IIPUT  "Enter  the  number  ot  the  selection  jrou  would  like.";  selX 
IF  selX  «  1  THEI 
CALL  DATARUI 
GOTO  100 
EID  IF 

IF  selX  «  2  THEI 
CALL  DATACOIV 

soon  2600,  10 

GOTO  100 
EID  IF 

IF  selX  «  3  THEI 
CALL  PLOTDATA 
GOTO  100 
EID  IF 
EID  SUB 


A-23 


SUB  MAXIMUM  (max,  DMA! !  ( ) ,  SI'/.) 
max  =  0 

FOR  X  =  1  TO  SIX 

DMAX(K)  =  CSIG(DMAX ! (X) ) 

NEXT  K 

FOR  K  =  1  TO  SIX 

IF  DMAX(K)  >  max  THEI  max  =  DMAX(K) 

HEXT  K 
END  SUB 

SUB  MINIMUM  (min,  DMI*!(),  SIX) 
min  =  0 

FOR  1  =  1  TO  SIX 

DMIN(X)  -  CSIG(DMII ! (X) ) 

NEXT  1 

FOR  1  =  1  TO  SIX 

IF  DMIN(l)  <  min  THEI  min  =  DMII(l) 

MEXT  1 
END  SUB 

SUB  PLOTDATA 

250  CALL  PLOTMEIU (TDCI ,  CHILI) 

CLS  0 
END  SUB 

SUB  PLOTMEIU  (TDCI ,  CHILI) 

DIM  TIME! (800),  PRES! (800),  P0SITI0IX(6) 

CLS  0 

VIEW  PRIIT  1  TO  25 
CHILI  =  "" 

PRIIT  PRIIT  "":  PRIIT  "" 

PRIIT  PRIIT  MM 

215  PRIIT  "Enter  test  mn  art  you  interested  in  or  hit  (Enter)  to" 
IIPUT  "return  to  the  Main  Mann";  TDCI 
IF  TDCI  »  ""  THEI 
CHILI  «  "7" 

GOTO  240 
END  IF 

200  CLS  0 

VIEW  PRIIT  1  TO  25 
TDCH  *  ""  +  TDCI  +  ".DAT" 

OPEI  "C:\SWS\DATAi\"  ♦  TDC1I  +  »»  FOR  IIPUT  AS  #5 
LIIB  IIPUT  15,  JUIXI 

IIPUT  *5,  POSITIOIX(I),  P0SITI0IX(2) ,  P0SITI0lX(3) , 
P0SITI0IX(4) ,  P0SITI0lX(6) ,  P0SITI0IX(6) 

CLOSE  *5 
PRIIT  "" 

PRIIT  "  Plot  Menu  Selections" 

PRIIT 

PRIIT  "1.  Plot  channel  1  (Position:  ";  POSITIOlX(l) ;  ")" 
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PRUT 

"2. 

Plot 

channel 

2 

(Position:  " 

P0SITI0iy.(2) 

II )  II 

PRUT 

"3. 

Plot 

channel 

3 

(Position:  " 

POSITIOIXO) 

M  ^  II 

PRUT 

"4. 

Plot 

channel 

4 

(Position:  " 

POSITIOIXO) 

II  J  M 

PRUT 

"5. 

Plot 

channel 

5 

(Position:  " 

POSITIOIXO) 

M  J  II 

PRUT 

"fl. 

Plot 

channel 

6 

(Position:  " 

POSITIOIXO) 

II  J  M 

PRUT 

"7. 

Enter  different  test  run." 

PRUT 

"8. 

Exit 

Plot  Program. " 

PRUT 

IIPUT  "Enter  the  item  number  you  would  like  to  view.  sel’/. 
IF  selX  =  1  THEM 

CHILI  =  LTRIM|(STR|(P0SITI0*X(1) )  ) 

IF  POSITIOIXO)  <  10  THE* 

CHILI  =  "0"  +  CHILI  ♦  »•■ 

EID  IF 
GOTO  225 
EID  IF 

IF  seiy.  =  2  THEI 

CHILI  =  LTRIHKSTRt  (P0SITI0I'/.(2)  )  ) 

if  posrrio*y.(2)  <  10  thei 

CHILI  =  "0"  +  CHILI  +  »" 

EID  IF 
GOTO  225 
EID  IF 

IF  aelX  =  3  THE* 

CHILI  =  LTRIMl(STRl(P0SITI0lX(3)  )  ) 

IF  POSITIOIXO)  <  10  THE* 

CHILI  a  "0"  +  CHILI  +  "" 

EID  IF 
GOTO  225 
EID  IF 

IF  eeiy.  =  4  THE* 

CHILI  =  LTRIHl(STRl(P0SITI0iy,(4)  )  ) 

IF  POSITIOIXO)  <  10  THE* 

CHILI  *  "0"  +  CHILI  ♦  "" 

EID  IF 
GOTO  225 
EID  IF 

IF  aelX  »  5  THEI 

CHILI  »  LTRIHI  (STRI  (POSITIOIXO) ) ) 

IF  POSITIOIXO)  <  10  THEI 
CHILI  »  "0"  +  CHILI  ♦  MM 
EID  IF 
GOTO  225 
EID  IF 

IF  aelX  >  6  THEI 

CHILI  «  LTRIHl ( STRl ( POS ITI OlX (6 ) ) ) 

IF  POSITIOIXO)  <  10  THEI 
CHILI  »  "0"  +  CHILI  ♦  "" 

EID  IF 
GOTO  225 
EID  IF 
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IF  selX  =  7  THEI 
GOTO  215 
EID  IF 

IF  selX  =  8  THEI 
CHILI  =  "7" 

GOTO  240 
EID  IF 

225  CALL  CHAIIEL(TDC|,  CHILI,  TIME! 0 ,  PRES!()) 

CALL  MAXIMUM (PMAX,  PRES!(),  TSTEPX) 

CALL  MI*IMUM(PMII,  PRES!(),  TSTEPX) 

CALL  MAXIMUM(TMAX,  TIME!(),  TSTEPX) 

CALL  MIIIMUMCTMII,  TIME!(),  TSTEPX) 

CLS  0 
SCREEI  9 

PRUT  TMAX ,  TMII ,  PMAX,  PMII 
AXIS  TMAX,  TMII,  PMAX,  PMII 
FOR  J  =  1  TO  TSTEPX 

IF  J  =  1  THEI  LIIE  (TIME! ( J) ,  PRES! (J))-(TIME! (J) ,  PRES!(J)) 
LIIE  -(TIME! (J) ,  PRES ! ( J) ) 

I  EXT  J 
00 

Al  =  IIKEYI 
LOOP  UITIL  Al  <>  "H 
GOTO  200 
240  EID  SUB 

SUB  VOLTAGE  (CH!(),  ARRAYXO ,  DIHAX,  IX,  FUGX.  TX,  TSX.  IlX,  I2X) 

I 

>  ************************************************************* 

'  *  SUBROUTIIE  VOLTAGE  * 

>  *  * 

'  *  This  subrout ins  takes  ths  binary  values  of  ths  sampled  * 

'  *  and  calculates  the  corresponding  voltages.  * 


FOR  I  =  t  TO  IX 
CH! (I)  *  0! 

I  EXT  I 

IF  «  »  1  THEI 
IlX  -  1 
12%  «  FUGX 
ELSE 

IlX  >  IlX  ♦  FLAGX 
12%  ■  12%  ♦  FUGX 
EID  IF 

FOR  J  «  1  TO  IX 

FOR  I  »  IlX  ♦  (J  -  1)  TO  I2X  +  (J  -  1)  STEP  1% 

CH! (J)  »  (2047  -  ARRAYX(I))  /  819.2  ♦  CH ! (J)  'COIVERT  TO 

VOLTAGES 

IEXT  I 
■EXT  J 
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PSIG  TRANSDUCER 


1000  FOR  I  *  1  TO  It 

CH!  (I)  =  CH !  ( I )  /  TSV. 
IEXT  I 

I 

END  SUB 
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Appendix  B.  Pressure  Transducer  Calibration  Curves 


The  calibration  curves  of  the  pressure  transducers  used  in  this  research  are 
presented  in  the  following  figures. 


I 

I 


VOLTS 


Figure  B.l.  Calibration  Curve  for  Pressure  Transducer  8510B-5  Serial  Number: 
PP81 


Figure  B.2.  Calibration  Curve  for  Pressure  Transducer  8510B-5  Serial  Number: 
79HB 
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Figure  B.3.  Calibration  Curve  for  Pressure  Transducer  8510B-5  Serial  Number: 
78HB 


10-0  —  TTTTTTTT-rn ~TTTTT-r  'I":  f  ;  I  I  r 1  'TT"’ T  m 

i  -i 

J  1 


Figure  B.4.  Calibration  Curve  for  Pressure  Transducer  8510B-15  Serial  Number: 
9BYP 
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VOLTS 

Figure  B.5.  Calibration  Curve  for  Pressure  Transducer  8510B-100  Serial  Number: 
35DF 


VOLTS 

Figure  B.6.  Calibration  Curve  for  Pressure  Transducer  8506B-5  Serial  Number: 
92BF 
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Figure  B.7.  Calibration  Curve  for  Pressure  Transducer  8506B-5  Serial  Number: 
68BF 


Figure  B.8.  Calibration  Curve  for  Pressure  Transducer  8506B-5  Serial  Number: 
74BF 


B-5 


Bibliography 


I 

I 

I 


1.  Lankford,  D.  W.  and  Rapaganani,  N.  L.  “Modeling  of  cavity  Flows  with 
Large  Base  Relief,”  International  Symposium  on  Gas-Flow  and  Chemical 
Lasers.  343-357.  Washington:  Hemisphere  Publications,  1979. 

2.  Korst,  H.,  H.  “A  Theory  for  Base  Pressures  in  Transonic  and  Supersonic 
Flow,”  Proceedings  of  the  National  Applied  Mechanics 

Division  Conference.  593-599.  New  York:  ASME  Press,  1956 

3.  Rudulf,  Hermann,  Dr.  Supersonic  Inlet  Diffusers  and  Introduction 

to  Internal  Aerodynamics.  Minneapolis:  Minneapolis-Honeywell  Regula¬ 
tor  Company,  1956. 

4.  Shapiro,  Ascher,  H.  The  Dynamics  and  Thermodynamics 

of  Compressible  Fluid  Flow.  New  York:  John  Wiley  anmd  Sons,  Inc., 
1953. 

5.  Hoffman,  Joe,  D.  and  Zucrow,  Maurice,  J.  Volume  I,  Gas  Dynamics.  New 
York:  John  Wiley  and  Sons,  Inc.,  1976. 

6.  Ames  Research  Staff.  Equations,  Tables,  and  Charts  for 
Compressible  Flow.  Report  No.  1135.  Ames  Aeronautical  Laboratory. 
Moffett  Field,  CA. 

7.  Zakanycz,  Stephen.  Turbulence  and  the  Mixing  of  Binary  Gases.  PhD 
dissertation.  Ohio  State  University,  OH,  1971. 

8.  Gates,  Capt  Gregory,  Engineer,  Research  Division.  Personal  interviews. 
HQ  AFWL,  Kirtland  AFB,  NM,  1  February  through  15  November  1989. 

9.  Air  Force  Weapons  Laboratory,  Air  Force  System  Command.  Contract 
F29601-80-C-0040  with  TRW  Corporation.  Drawing  X430038,  Sheet  2. 

10.  Qua  Tech.  Users  Manual  for  the  Multichannel  Simultaneous  Sampling  Sys¬ 
tem.  Qua  Tech,  Inc,  Akron,  OH,  undated. 

11.  ENDEVCO.  General  Catalog.  ENDEVCO  San  Juan  Capistrano,  CA, 
1988. 

12.  AMETEK  Mansfield  and  Green  Division.  Handbook  Operating 

And  Service  Instructions.  Form  No.  75-109  (Rev.  4).  AMETEK  Mansfield 
and  Green  Division,  Largo,  FL,  June  1984. 

13.  Holder,  D.,  W.  and  North,  R.,  J.  AGARDograph  23,  Optical  Methods  for 
Examining  the  flow  in  High-speed  Wind  Tunnel,  Part  I  Schlieren 
Methods.  North  Atlantic  Treaty  Organization,  Avisory  Group  For  Aero¬ 
nautical  Research  and  Development.  1956. 


BIB-1 


14.  Tanis,  Capt  Frederick,  PhD  Canidate.  Personal  interviews,  School  of  En¬ 
gineering,  Air  Force  Institute  of  Technology,  Wright-Patterson  AFB  OH, 
4  January  through  15  November  1989. 


Vita 


Captain  Stephen  W.  Stiglich, 

Pennsylvania.  At  the  age  of  fourteen  he  accepted  Christ. 


n  Au¬ 


gust  1971  he  enlisted  in  the  United  States  Air  Force.  He  received  an  instrument 
rated  commercial  pilots  license  in  January  1976.  .  He  was  selected  for  the  Air¬ 


man  Education  and  Commissioning  Program  (AECP  )  in  October  1980,  and 
he  entered  the  University  of  Texas  at  Austin  in  May  1981.  He  graduated  from 
the  University  of  Texas  at  Austin  with  a  Bachelor  of  Aerospace  Engineering 
degree  in  December  1984.  Following  graduation  he  attended  Officer  Training 
School  and  was  commissioned  in  April  1985.  His  first  duty  station  was  the 
Defense  Dissemination  Program  Office  at  Space  Division  where  he  served  as 
a  Special  Projects  Engineer  for  three  years.  Captain  Stiglich  entered  the  Air 


Force  Institute  of  Technology  in  May  1988. 


VITA-1 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 

AFTT/GAE/ENY/89D-36 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


Form  Approved 
OMB  No.  0704-0188 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

School  of  Engineering 


6b.  OFFICE  SYMBOL 
(If  applicable) 

AFIT/ENY 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Air  Force  Institute  of  Technology  (AU) 
Wricrht-Patterson  AFB  OH  45433-6583 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


8a.  NAME  OF  FUNDING  /SPONSORING 

8b.  OFFICE  SYMBOL 

ORGANIZATION 

(If  applicable) 

Air  Force  Weapons  Laboratory 

AFWL/ARDK 

Sc.  ADDRESS  (City,  State,  and  ZIP  Code) 

Air  Force  Weapons  Laboratory 
Kirtland  AFB  NM  87117 

10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 


PROJECT 

TASK 

NO. 

NO 

1 1 .  TITLE  (Include  Security  Classification) 

EXPERIMENTAL  INVESTIGATION  OF  A  CHEMICAL  LASER  CAVITY  FLOWFIEXD 


WORK  UNIT 


12.  PERSONAL  AUTHOR(S) 

Stephen  W.  Stiglich,  Jr. ,  B.S. ,  Capt,  USAF 


13a.  TYPE  OF  REPORT  13b  TIME  COVERED 

MS  Thesis  from _ to _ 


16.  SUPPLEMENTARY  NOTATION 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  IIS.  PAGE  COUNT 

1989  December  I  120 


FIELD 

GROUP 

20 

04 

COSATI  CODES 


SUB-GROUP 


18.  SUBJECT  TERMS  { Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Laser  Cavities  Chemical  Lasers  Gas  Dynamics 
Flow  Fields  Base  Flew  Flew  Visualization 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Thesis  Advisor;  William  C.  Elrod 
Professor 

Department  of  Aeronautics  and  Astronautics 


20.  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT. 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

William  C.  Elrod,  Professor 


DD  Form  1473,  JUN  86  a 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
□  DTIC  USERS  UNCLASSIFIED  _ 


22b.  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

513)  255-3517  ENY 


III 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


Abstract 

Chemical  lasers  require  a  cavity  that  establishes  and  maintains  the  proper  gas 
dynamic  properties  during  lasing.  The  design  and  performance  of  a  flow  system 
capable  of  supporting  the  hypersonic  flow  conditions  in  a  lasing  cavity  are  described. 
Using  cold  air  as  the  working  medium,  the  flow  control  system  configuration  and 
nozzle-cavity-supersonic  diffuser  assembly  configuration  were  developed  to  establish 
acceptable  flow  conditions  in  the  test  section.  Performance  evaluation  was  based 
on  pressure  measurements  in  the  nozzle-cavity-diffuser  assembly  and  schlieren  pho¬ 
tographs  of  the  flowfield  in  the  cavity. 

Flow  conditions  in  the  test  section  were  broken  up  into  three  different  regions: 
flow  in  the  hypersonic  nozzles,  flow  in  the  base  region  and  flow  in  the  cavity  region. 

Flow  in  the  nozzles  was  analyized  using  one-dimensional,  steady,  isentropic 
flow  theory.  Test  results  indicated  that  the  hypersonic  nozzles  performed  to  design 
specifications. 

The  Korst  two-dimensional  base-pressure  flow  model  was  used  to  describe 
the  flow  in  the  nozzle  exit  plane  and  base  region.  Experimentally  calculated  Mach 
numbers  and  static  pressures  corresponded  very  closely  to  theoretical  values. 

Static  pressure  ports  and  schlieren  photographs  were  used  to  describe  the  flow- 
field  conditions  in  the  cavity  region.  Pressure  measurements  indicated  that  super¬ 
sonic  conditions  were  reached  in  the  cavity  for  specific  supersonic  diffuser  throat 
areas  settings,  but  conditions  were  short  lived.  Boundary  layer,  frictional,  and 
three-dimensional  effects  were  suspected  as  the  main  contributors  to  the  flowfield 
degradation. 


UNCLASSIFIED 


